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Rational  Techniques  for  Evaluating  Potential  Sands  for  Beach  Nourishment  (CR 
DRP-93-2) 


ISSUE:  The  rational  design  of  beach- 
nourishment  projects  requires  the  ability  to 
calculate  the  final  configuration  of  the  added 
sand  volume.  This  capability  is  essential  for 
quantitative  evaluation  of  the  relative  merits 
of  various  borrow  areas  and  in  benefit/cost 
analysis  of  such  projects. 

RESEARCH:  A  methodology  was  devel¬ 
oped  for  predicting  the  equilibrium  beach  pro¬ 
file  resulting  from  placemen:  of  an  arbitrary 
volume  of  material  with  an  arbitrary  grain- 
size  distribution  on  a  profile  of  arbitrary  shape 
and  grain-size  distribution.  Tf'e  methodology 
depends  on  the  theory  of  equilibrium  profile 
shape  and  is  proposed  as  an  alternative  to  tra¬ 
ditional  compatibility  and  overfill  ratio  factors 
for  borrow  and  native  material. 

SUMMARY:  Various  available  methods  for 
relating  the  overall  qualities  of  borrow  and 
native  sediments  were  reviewed.  These  meth¬ 
ods  generally  focus  on  comparing  grain-size 
characteristics  rather  than  response  to  a  wave 
and  tide  regime.  Thus  the  procedures  must  be 
considered  ad  hoc  and  the  results  not  truly  rep¬ 
resentative  of  a  beach  nourishment  material. 
Also,  the  methods  cannot  be  used  to  deter¬ 


mine  additional  dry  beach  width,  a  factor  nec¬ 
essary  in  benefit/cost  analysis  and  for  project 
design  and  planning  operations. 

A  methodology  and  a  computer  program  were 
developed  for  predicting  beach  shapes  rele¬ 
vant  to  beach  nourishment  using  sediments  of 
arbitrary  sorting.  The  theory  was  then  applied 
in  four  specific  examples  covering  a  range  of 
beach,  native  sediment,  and  fill  material  condi¬ 
tions.  Two  sets  of  field  data  were  located  that 
provided  comparison  with  and  evaluation  of 
the  general  methodology. 

The  method  presented  in  the  report  is  consid¬ 
ered  by  the  authors  to  be  one  step  toward  a  ra¬ 
tional  procedure  for  assessing  the  complex 
performance  of  nourishment  projects  with  re¬ 
alistic  sediment  characteristics. 

AVAILABILITY  OF  REPORT:  The  report 
is  available  through  the  Interlibrary  Loan  Ser¬ 
vice  from  the  U.S.  Army  Engineer  Waterways 
Experiment  Station  (WES)  Library,  telephone 
number  (601 )  634-2355.  National  Technical 
Information  Service  (NTIS)  report  numbers 
may  be  requested  from  WES  Librarians.  To 
purchase  a  copy  of  the  report,  call  NTIS  at 
(703)  487-4780. 
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Summary 


Rational  design  of  beach  nourishment  projects  requires  the  ability  to 
calculate  the  geometry  of  the  added  sand  volume.  This  capab;iity  is  essential 
for  quantitative  evaluation  of  the  relative  merits  of  various  borrow  areas  and 
for  benefit/cost  analysis  of  such  projects.  In  many  cases,  the  material  may  be 
a  by-product  of  a  dredging  project  carried  out  for  purposes  other  than  beach 
nourishment,  and  the  dredged  material  to  be  placed  may  have  a  different  grain 
size  distribution  than  the  original  (native)  beach. 

This  report  presents  a  new  methodology  for  predicting  the  equilibrium 
beach  profile  resulting  from  placement  of  an  arbitrary  volume  of  material  with 
an  arbitrary  grain  size  distribution  on  a  profile  of  arbitrary  shape.  The 
methodology  developed,  which  depends  on  the  theory  of  equilibrium  profile 
shape,  is  proposed  as  an  alternative  to  traditional  compatibility  and  overfill 
ratio  factors  for  borrow  and  native  material.  The  methodology  considers  two- 
dimensional  (cross-shore)  conditions. 

The  theory  is  first  developed  and  characteristics  of  equilibrium  beach 
profiles  relevant  to  beach  nourishment  projects  are  presented.  The  theory  is 
then  applied  in  four  specific  examples  covering  a  range  cf  beach,  native 
sediment,  and  fill  material  conditions.  A  computer  program  developed  for 
applying  the  method  is  given  in  an  appendix. 
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1  Introduction 


Rational  design  of  beach  nourishment  projects  requires  the  ability  to 
calculate  the  geometry  of  the  added  sand  volume,  both  as  a  function  of  space 
and  time.  This  capability  is  essential  for  quantitative  evaluation  of  the 
relative  merits  of  various  borrow  areas  and  in  benefit/cost  analysis  of  such 
projects,  including  the  volumes  and  timing  of  renourishments.  In  many  cases 
the  material  may  be  a  by-product  of  a  dredging  project  carried  out  for  other 
purposes.  Particular  design  aspects  of  significance  include  the  equilibrated 
beach  profile,  especially  the  additional  dry  beach  width  and  the  longevity  of 
the  project.  Traditionally,  attempts  to  quantify  benefits  of  beach  nourishment 
have  utilized  "compatibility"  and  "overfill"  factors,  which  relate  the  sand  sizes 
of  borrow  and  native  material. 

In  a  general  sense,  the  problem  of  the  evolution  of  beach  nourishment 
projects  can  be  considered  as  occurring  over  two  more  or  less  distinct  time 
scales  (see  Figure  1).  The  beach  profile,  which  is  usually  placed  at  a 
relatively  steep  slope,  equilibrates  over  a  fairly  short  time  scale,  perhaps  with 
a  "folding  time"  on  the  order  of  several  years.  The  time  scale  associated  with 
the  planform  evolution  depends  primarily  on  the  length  of  the  project  and  the 
wave  climate;  for  longer  projects  (say  greater  than  several  kilometers),  the 
time  scale  is  on  the  order  of  decades.  Regardless  of  whether  these  two  time 
scales  are  distinct  or  not,  it  is  useful  to  consider  the  equilibrium  beach  profile 
associated  with  the  volume  and  texture  (i.e.,  grain  size  distribution)  of  the 
nourishment  material. 

This  report  presents,  for  two-dimensional  conditions,  procedures  for 
predicting  the  equilibrium  beach  profile  resulting  from  placement  of  an 
arbitrary  volume  of  material  with  an  arbitrary  grain  size  distribution  on  a 
profile  of  arbitrary  shape. 

This  report  is  organized  as  follows:  Chapter  2  reviews  the  background 
relative  to  efforts  to  quantify  the  suitability  and/or  effectiveness  of  materials 
for  beach  nourishment.  Additionally,  characteristics  of  equilibrium  beach 
profiles  relevant  to  beach  nourishment  are  presented.  These  latter  results 
pertain  to  perfectly  sorted  sediments.  Chapter  3  describes  and  illustrates  with 
examples  the  methodology  developed  in  the  present  study  for  calculating 
equilibrium  profiles  for  sediments  of  arbitrary  sorting.  A  computer  program 
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Original  Shoreline 


'Spreading  Out"  Losses 


Sand  Moves  Offshore  to 
Equilibrate  Profile 


Nourished  Shoreline 


'Spreading  Out"  Losses 


a)  Plan  View  Showing  "Spreading  Out"  Losses 
and  Sand  Moving  Offshore  to  Equilibrate  Profile 


Dry  Beach  Width 
(Fine  Sand) 


Dry  Beach  Width  (Coarse  Sand) 
initial  Placed  Profile 

Equilibrated  Profile  (Coarse  Sand) 
Sea  Level 


Original  Profile 


Equilibrated  Profile  (Fine  Sand) 


b)  Elevation  View  Showing  Original  Profile,  Initial  Placed  Profile, 
and  Adjusted  Profiles  That  Would  Result  from  Nourishment 
Project  with  Coarse  and  Fine  Sands 

Figure  1 .  Sand  transport  losses  and  beach  profiles  associated  with  a  beach 
nourishment  project 
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developed  for  this  purpose  is  described  in  terms  of  the  general  algorithms 
employed.  Chapter  4  presents  the  results  of  a  limited  set  of  laboratory  studies 
carried  out  in  conjunction  with  this  study.  Chapter  5  describes  relevant  field 
data  from  Delray  Beach,  Florida,  and  Jupiter  Island,  Florida.  Both  areas  have 
experienced  multiple  beach  nourishment  programs.  Chapter  6  provides  the 
summary,  conclusions,  and  recommendations  for  further  research.  Appendix 
A  is  a  list  of  the  computer  program  provided  to  carry  out  calculations  of 
equilibrium  profiles  for  sediments  of  arbitrary  sorting.  Additionally,  listings 
of  the  input  and  output  files  for  two  of  the  examples  presented  in  this  report 
are  provided.  Appendix  B  provides  a  detailed  discussion  of  the  computer 
program.  Appendix  C  contains  additional  data  for  the  Delray  Beach,  Florida, 
nourishment  project.  Appendix  D  is  a  notation  of  symbols  used  in  this  report. 
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2  Background 


Various  investigators  have  proposed  procedures  for  relating  the  overall 
qualities  of  borrow  and  native  sediments.  As  a  rule  of  thumb,  an  attempt  is 
made  to  locate  borrow  material  with  granulometric  characteristics  similar  to 
those  of  the  native  material.  Generally,  the  available  procedures  focus  on 
comparing  the  grain  size  characteristics  rather  than  on  their  response  to  a 
wave  and  tide  regime.  Thus  these  methods  must  be  considered  as  ad  hoc  and 
not  truly  representative  of  the  performance  of  a  beach  nourishment  material. 
Moreover,  these  procedures  do  not  address  the  additional  dry  beach  width,  a 
factor  of  primary  concern  to  the  designer  and  funding  entities.  A  short  review 
of  the  various  available  methods  follows. 

Krumbein  and  James  (1965)  proposed  a  method  which  considered  the  grain 
size  distributions of  the  borrow  and  native  materials  to  be  each  repre¬ 
sented  by  log-normal  distributions  as  proposed  earlier  by  Krumbein  (1957) 


v/^o 

in  which  <l>  is  the  sediment  diameter  expressed  in  phi  units,  defined  as 


<ti=-logj(J)(mm)) 


(2) 


where,  as  indicated,  D  is  the  sediment  diameter  in  millimeters  and  jx  and  a  are 
the  sample  mean  and  standard  deviation  in  phi  units.  This  method  defined 
compatibility  of  the  borrow  material  on  the  basis  of  the  proportion  of  borrow 
material  distribution  that  was  common  with  the  native  sand  size  distribution. 
This  approach  appears  somewhat  reasonable  in  discounting  the  finer  fraction 
of  the  borrow  material,  but  less  reasonable  in  discounting  similarly  the  pro¬ 
portion  of  coarse  material  that  is  in  excess  relative  to  the  native  sand. 


’  For  convenience,  symbols  and  abbreviations  are  listed  in  the  notation  (Appendix  D). 
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James  (1974)  developed  a  complex  method  addressing  the  relative 
renourishment  frequency  for  different  sand  characteristics;  however,  this 
procedure  only  considered  longshore  sediment  transport  and  considered  the 
nourishment  project  to  be  located  in  an  area  where  the  ambient  longshore 
sediment  transport  had  been  interrupted  completely. 

Dean  (1974)  presented  a  method  that  attempted  to  address  the  deficiency 
(noted  above)  of  the  earlier  Krumbein  and  James  method.  The  borrow 
material  was  only  discounted  for  the  excessive  proportion  of  fines  present; 
excessive  proportions  of  coarser  material  were  included  in  the  compatible 
fraction.  However,  it  was  considered  that  all  the  fine  fraction  smaller  than  a 
critical  value  was  lost.  This  method  resulted  in  a  considerably  higher 
compatibility  than  that  of  Krumbein  and  James  (196S). 

James  (1975)  developed  a  renourishment  factor  based  on  the  relative  char¬ 
acteristics  of  the  borrow  and  native  sand  characteristics.  Similar  to  earlier 
methods,  this  procedure  was  based  on  size  distributions  rather  than  associated 
equilibrium  profiles.  Compared  to  Dean’s  method  (1974),  the  primary 
difference  is  the  retention  of  a  portion  of  the  fine  fraction  in  the  compatibility 
considerations.  At  present,  the  methods  of  James  (1975)  are  those 
recommended  in  the  Shore  Protection  Manual  (1984). 

Companion  to  the  problem  of  defining  borrow  material  compatibility  is  that 
of  sampling  across  the  pre-nourishment  profile  to  establish  the  "native”  sand 
characteristics.  This  problem  has  been  addressed  by  several  investigators. 

Although  the  primary  focus  of  the  present  report  is  the  equilibrium  beach 
profile  for  sand  of  arbitrary  distribution,  it  is  useful  to  consider,  for  back¬ 
ground  purposes,  profiles  which  result  for  the  idealized  case  of  uniform 
borrow  and  native  sediment  sizes. 

Dean  (1991)  has  considered  equilibrium  beach  profiles  represented  by 
h(y)  =  first  proposed  by  Bruun  (1954)  and  found  later  by  Dean  (1977)  in 
an  analysis  of  more  than  500  profiles  extending  from  the  eastern  end  of  Long 
Island  around  Florida  to  the  Gulf  of  Mexico  border.  In  this  equation  h  is 
water  depth,  y  is  distance  seaward  from  the  shoreline,  and  A  is  an  empirical 
coefficient  called  the  scale  parameter  or  simply  "A  value."  Moore  (1982) 
investigated  the  relationship  between  the  sediment  scale  parameter  A  and  the 
sand  diameter  D,  and  established  the  results  shown  by  the  curved  line  in 
Figure  2.  Later  Dean  (1987)  simply  transformed  this  A  versus  D  relationship 
to  A  versus  w  where  w  is  the  fall  velocity.  The  result  was  found  to  be  well- 
approximated  by  the  straight  line  in  Figure  2. 

It  has  been  shown  that  three  types  of  profiles  could  occur,  depending  on 
the  relative  sizes  of  the  borrow  and  native  sands.  These  are  termed  "intersect¬ 
ing,"  "non-intersecting,"  and  "submerged"  profiles  and  are  illustrated  in 
Figure  3.  The  reader  is  referred  to  the  paper  of  Dean  (1991)  for  the  criteria 
separating  the  three  profile  types  and  the  volumes  required  to  achieve,  for 
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example,  a  desired  additional  dry  beach  width  of  the  nourished  profile  (for 
intersecting  and  non-intersecting  profiles). 

The  significance  of  these  three  profile  types  can  be  seen  by  referring  to 
Figures  4,  5,  and  6.  Figure  4  shows  the  effect  on  additional  dry  beach  width 
of  placing  the  same  volume  (340  m^/m)  of  sand  of  four  different  grain  sizes. 
In  the  upper  panel,  the  sediment  is  coarser  than  the  native  and  the  profiles 
intersect  with  an  additional  dry  beach  width  of  92.4  m.  Panel  b  shows  the 
effect  of  using  sand  of  the  same  size  as  the  native  resulting  in  a  dry  beach 
width  of  45.3  m.  Panels  c  and  d  present  the  results  for  decreasing  sediment 
size;  in  Panel  d,  the  dry  beach  width  is  zero. 

Figure  5  shows  the  effect  of  nourishing  with  various  quantities  of  a 
sediment  that  is  smaller  than  the  native.  With  the  same  sediment  size,  the 
volumes  of  sediment  increase  from  Panels  a  to  d.  With  increasing  volume, 
the  landward  and  seaward  extent  of  the  nourished  regions  increase  and  in 
Panel  d,  sufficient  volume  has  been  added  to  achieve  a  transition  between 
submerged  and  non-intersecting  profiles. 

These  types  of  results  can  also  be  presented  as  shown  in  Figure  6  for  an 
example  that  is  in  a  form  more  representative  of  beach  nourishment  concerns. 
This  figure  shows  the  relationship  of  additional  dry  beach  width  Ay  versus 
volume  added  V  for  three  values  of  nourishment  sediment  scale  parameter 
Other  variables  common  to  the  three  cases  are:  berm  height,  B  =  2.0  m, 
depth  of  effective  motion,  h,  =  8.0  m,  and  native  sediment  scale  parameter, 
Af,  =  0.1  m‘'^  Of  interest  is  that  for  A,:  =  0  12  >  A^,  the  profiles  are 
initially  intersecting  and  the  additional  dry  beach  width  increases  relatively 
rapidly.  However,  with  increasing  volume,  the  profile  becomes 
non-intersecting  and  the  slope  d(,Ay)ldV  is  approximately  the  same  as  that  for 
Af  =  A^,  which  is  almost  a  constant.  For  A,,  =  0.08  m"^  <  Af^,  for  small 
volumes  of  sediment,  there  is  no  additional  beach  width,  i.e.,  the  profile  is 
submerged.  However,  with  increasing  volumes,  there  is  a  critical  volume  at 
which  the  landward  end  of  the  submerged  profile  just  reaches  the  shoreline; 
for  still  greater  volumes,  the  profile  becomes  a  non-intersecting  profile  and 
remains  so  for  increasing  volumes.  For  this  case,  the  slope  d{Ay)ldV  is 
essentially  constant  and  parallel  to  the  case  of  Af^  =  It  is  stressed  that  all 
of  these  results  apply  for  equilibrium  profiles  and  that  the  equilibration 
process  may  take  several  years  to  complete.  In  most  nourishment  projects, 
the  sand  is  placed  steeper  than  equilibrium  and  will  provide  greater  additional 
dry  beach  width  than  equilibrium  during  the  equilibration  process. 

In  general  it  can  be  shown  (Dean  1991)  that  the  non-dimensional  additional 
beach  width  AyAV.  is  related  to  the  non-dimensional  volume  added  K/flW., 
non-dimensional  berm  height  BIh,,  and  ratio  of  fill  to  nourishment  sediment 
scale  factors  AflAf^,  i.e.. 


6 


Chapter  2  Background 


(3) 


Az  -  f(  ^  A 

W,  ’  A^j 

where  h,  and  W,  are  the  depth  of  limiting  motion  and  width  of  the  active 
profile  for  the  initial  native  profile.  Figures  7  and  8  present  this  relationship 
for  B/h,  =  1/2  and  1/4,  respectively.  Several  features  of  these  figures  are  of 
interest.  First,  for  >  1.2  (approximately),  there  is  little  additional  dry 
beach  width  gained  for  coarser  (greater  A^JA,,)  material  used.  Secondly,  for 
AfJA^  <  1 ,  and  a  fixed  volume,  there  is  a  rapidly  decreasing  dry  beach  width 
with  decreasing  AflA,/.  Finally,  the  transition  from  intersecting  to  non¬ 
intersecting  profiles  is  indicated  by  the  bold  line  in  Figures  7  and  8  and  the 
transition  from  non-intersecting  profiles  to  submerged  profiles  occurs  at  the 
vertical  asymptotic  lines  at  the  left  end  of  each  of  the  curves. 
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SEDIMENT  FALL  VELOCITY,  w  (cm/s) 


ariation  of  sediment  scale  parameter  A  with  sediment  size  and  fall  velocity 


Figure  3.  Three  generic  types  of  nourished  profiles:  (a)  intersecting  profile, 
(b)  non-intersecting  profile,  and  (c)  submerged  profile 


Chapter  2  BackgrourKi 


9 


ELEVATION 


Intersecting  Profiles, 

0.1m’ ^^.Ap  =  0.14m’ 
45.3m 
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Figure  4.  Four  examples  of  decreasing  Af  with  same  added  volume  per  unit 
beach  length 
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INCREASING  VOLUMES  OF  ADDED  SAND  ELEVATION  (m) 
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Figure  5.  Effect  of  increasing  volume  of  sand  added  on  resulting  beach  profile 
>4^=0. 1  m’'^,  -4^  =  0. 2  m’'®,  h.  =  6.0  m,  8=1.5  m 
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A'  =  Ap/A^ 


Figure  7.  Variation  of  non-dimensional  shoreline  advancement  Ay/UV.  with  A’ 
and  V.  Results  shown  for  h.lB  =  2.0 
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A'  =  Af/An 


Figure  8.  Variation  of  non-dimensional  shoreline  advancement  LytW.  with  A' 
and  V.  Results  shown  for  h.lB  =  4.0 
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3  Methodology  Based  on 
Equilibrium  Beach  Profiles 


The  method  developed  herein  is  based  on  an  equilibrium  profile  for  a 
sediment  volume  with  a  distribution  of  sediment  sizes.  The  following 
assumptions/considerations  are  made: 

a.  A  volume  of  sand  V  per  unit  beach  length  is  placed  at  a  slope  steeper 
than  equilibrium  in  conjunction  with  a  beach  nourishment  project. 

b.  The  sand  is  well-mixed  at  the  time  of  placement. 

c.  This  sand  will  be  reworked  such  that  the  volume  removed  from  the 
placement  cross  section  is  sufficient  to  extend  the  nourished  equilibrium 
profile  out  to  a  specified  depth  A.  of  limiting  motion  (see  Figure  9)  or 
to  its  intersection  with  the  initial  profile. 

d.  Within  the  zone  of  sediment  removal  from  the  placement  cross  section, 
sorting  occurs  down  to  a  specified  thickness 

e.  The  available  sediment  is  sorted  across  the  profile  with  the  coarser 
fraction  remaining  in  the  berm  and  shallower  water  region  and 
the  finer  sediment  distributed  offshore. 

/.  The  volumes  of  sediment  removed  and  deposited  are  equal;  i.e.,  volume 
is  conserved. 

With  the  above  basis,  the  procedure  can  be  considered  as  one  of  locally 
establishing  segments  of  an  equilibrium  profile  consistent  with  the  local  A 
value  and  of  balancing  sediment  volumes.  Because  the  equilibrium  beach 
profile  form  h  =  Ay^'^  yields  an  unrealistic  infinite  slope  at  y  =  0,  the 
modified  form  was  used,  which  recognizes  the  effect  of  gravity  for  the  larger 
slopes 
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(4) 


5 


as  initially  proposed  by  Dean  (1983)  and  later  shown  by  Larson  (1988)  and 
Larson  and  Kraus  (1989)  to  be  derivable  from  the  breaking  wave  model  of 
Dally,  Dean,  and  Dalrymple  (1985)  under  the  consideration  of  uniform  wave 
energy  dissipation  per  unit  volume.  In  Equation  4,  s  is  the  beach  face  slope. 
It  can  be  shown  easily  that  in  shallow  water. 


i.e.,  the  beach  is  planar,  consistent  with  measurements  in  nature.  In  deeper 
water.  Equation  4  approximates  h  = 

Because  the  A  value  is  now  local,  the  depth  at  a  location  y  +  dy  is 
referenced  to  the  depth  at  y  based  on  Equation  4, 

hiy^dy)  =  h(y)  +  (6) 

(dyidh) 


where 


s  2\a^] 

/*(y) 


A  <  0 
h>0 


(7) 


and  the  dy  values  are  maintained  reasonably  small,  on  the  order  of  1-2  m. 

A  step-by-step  discussion  of  the  procedure  is  as  follows  and  is  illustrated  in 
the  program  flowchart,  Figure  10. 


a.  With  the  specified  initial  profile  /io(y),  added  volume  V,  berm  height  B, 
and  placemen?  slope  s^,,  the  placed  profile  h^fy)  is  determined  by 
iteration  such  'J-ni  the  volume  out  to  the  location  where  h^iy)  =  ^(y)  is 
the  volume  piu  ’.ed.  This  procedure  also  determines  the  berm  advance¬ 
ment  Ay„. 

b.  Trial  values  of  the  volume  sorted  and  equilibrated  berm  advance¬ 
ment  DYEQ  are  assumed  (refer  to  Figure  9  for  definition  of  variables). 
For  each  pair  of  these  quantities,  the  equilibrium  profile  is  advanced 
from  y  to  y  -t-  dy,  where  dy  is  constant,  say  1-2  m.  This  advancement 
is  in  accord  with  Equations  6  and  7.  The  local  A  is  that  associated  with 
the  diameter  for  the  coarser  fraction  of  the  sediments  that  has  not  been 
deposited  up  to  y  in  the  equilibration  process  (see  Figure  1 1).  This 
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Figure  10.  Flow  diagram  for  problem  solution 
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Vqen  =  Horizontal  Inclined  Hatched  Volumes 
VusED  =  Vertical  +  Inclined  Hatched  Voumes 


NVHl  dSIlVHS  NOIlUOdOUd 
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Figure  1 1 .  Method  of  determining  grain  size  for  computation  of  current  increment  of  equilibrium  profile 


step-by-step  advancement  is  continued  until  the  depth  equals  the 
specified  terminal  depth  h,  or  until  the  equilibrium  profile  intersects 
with  the  initial  profile.  At  that  stage,  the  volume  actually  generated 
through  erosion  of  the  placed  profile  is  substituted  for  the  volume 
available  and  the  equilibrated  berm  advancement  DYEQ  is  held  fixed  in 
this  inner  loop  and  the  process  repeated.  This  inner  loop  (with  DYEQ) 
fixed  is  repeated  until  values  in  two  successive  iterations  agree 
within  an  acceptable  limit. 

c.  The  value  of  DYEQ  is  changed  to  attempt  to  ensure  that  the  associated 
value  of  h  =  h.  or  profile  intersection  will  be  achieved  coincident  with 
the  deposition  of  VGEN  for  that  value  of  DYEQ.  The  DYEQ  at  the 
/c-t- 1  iteration  is  based  on  the  following  simple  algorithm 


DY£<?**‘  =  DYEQ^  +  f*"‘(ADY£Q) 


(8) 


in  which  ADYEQ  is  specified  as  some  reasonable  value,  say  2  or  5  m, 
and  F'  =  1,  £*  =  4il  for  k=2  and  the  positive  and  negative  signs 

apply  depending  on  whether  Voen  >  (£^  =  1)  or 

(F^  =  -1).  In  subsequent  iterations  (it>2),  £**'  =  £*  if  the  sign  of 
'^oEN  ■  '^usED  dW  not  change  in  the  preceding  iteration  and 
=  ^.5  f*  if  a  sign  change  did  occur. 

It  is  noted  that  the  solution  procedure  structure  is  identical  for  both  ideal¬ 
ized  and  arbitrary  grain  size  distributions.  In  this  context,  idealized  refers  to 
grain  size  distributions  given  by  Equation  1.  Additionally,  the  method  can  be 
applied  for  arbitrary  initial  profiles. 

Examples 


Methods  developed  in  the  earlier  sections  of  this  paper  will  be  illustrated 
with  examples. 

Example  1 

Idealized  initial  profile  and  log-normal  size  distribution,  non-intersecting 
profile; 

In  this  example,  the  initial  profile  was  specified  as  characterized  by  the 
following: 

Uniform  sand  size:  A  =  0. 1  m"-’  (D  =  0.20  mm) 

Berm  height;  B  =  1.5  m 
Beach  face  slope;  s„  =  1:10 
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The  characteristics  of  the  nourishment  material  are  as  follows: 

Volume  added:  140.0  mVm 

Lx)g-normal  sand  size:  fi  =  1.6<^  (D  =  0.33  mm),  a  =  0.40<^ 

Berm  height:  B  =  1.5  m 

Placed  slope  ;  =  1:10 

Equilibrium  beach  face  slope:  Jeq  =  l-‘20 

Mixed  depth:  =  0.2  m 

Depth  of  active  motion:  A.  =  6  m 

Figure  12  presents  the  size  distribution  of  the  sediment.  The  volume  and 
slope  above  yielded  a  placed  shoreline  advancement  Ay„  of  S8.2  m. 

Figure  13  shows  the  initial,  placed,  and  equilibrium  beach  profiles.  It  is 
seen  that  for  this  case  the  equilibrated  shoreline  advancement  is  20.3  m.  The 
volume  eroded  for  this  case  is  66.8  m^  and,  of  course,  the  volume  eroded  is 
equal  to  the  volume  deposited.  For  this  example,  the  equilibrium  profile 
extends  to  an  offshore  distance  of  SSO  m,  where  it  reaches  the  specified  depth 
of  6  m. 

Example  2 

Idealized  initial  profile  and  log-normal  size  distribution,  intersecting 
profiles: 

The  characteristics  of  this  example  are  the  same  as  for  Example  1  except 
the  sorting  coefficient  a  of  the  placed  sand  is  0.1<^  rather  than  0A<t>.  The 
cumulative  sediment  size  distribution  is  shown  in  Figure  12.  For  this  case  the 
equilibrium  and  original  profiles  intersect  at  a  depth  of  4.69  m,  which  is 
located  at  a  distance  of  425  m  offshore.  Hie  initial,  placed,  and  equilibrium 
profiles  are  presented  in  Figure  14. 

Example  3 

Sand  smaller  than  native,  near-zero  shoreline  advancement: 

The  initial  profile  was  specified  as  follows: 

Uniform  sand  size:  (D  =  0.2  mm) 

Berm  height:  B  =  1.5  m 
Beach  face  slope:  s„  =  1:10 

The  characteristics  of  the  nourished  profile  are: 

Volume  added:  240  m^/m 

Log-normal  sand  size:  fi  =  2.64<^  (D  =  0.16  mm),  o  =  0. 10<^  (see 
Figure  12) 

Berm  height:  B  =  1.5  m 
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Figure  12.  Cumulative  grain  size  distribution,  Examples  1  through  4 
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Figure  13.  Example  1,  original,  placed,  and  equilibrium  profiles.  Case  of  non-intersecting  profiles.  Idealized  grain  size 
distribution 
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Figure  14.  Example  2,  original,  placed,  and  equilibrium  profiles.  Case  of  intersecting  profiles.  Idealized  grain  size  distribution 


Placed  slope:  s,  =  1:10 
Equilibrium  beach  face  slope:  igg  =  1:10 
Mixed  depth:  =  0.2  m 

Depth  of  active  motion:  A.  =  6  m 

This  example,  presented  in  Figure  15,  illustrates  conditions  near  a  transition 
to  a  submerged  profile. 

Example  4 

User-specified  initial  beach  profile  and  sediment  characteristics: 

In  this  example,  the  initial  beach  profile  was  specified  by  nine  points. 
Between  these  points,  the  profile  is  considered  as  a  series  of  straight  line 
segments.  The  sediment  size  distribution  is  specified  as  linear  as  shown  in 
Figure  12.  A  sediment  volume  of  600  m?lm  has  been  added.  Other  variables 
are  similar  to  those  specified  in  Example  3  and  are  shown  along  with  the 
initial,  placed,  and  equilibrium  profiles  in  Figure  16.  The  equilibrium  profile 
is  of  the  non-intersecting  type;  however,  intersection  nearly  occurs.  Other 
characteristics  are  similar  to  those  in  the  three  previous  examples. 

Example  5 

Variation  of  additional  dry  beach  width  versus  volume  added: 

Figure  17  presents  the  variation  of  additional  dry  beach  width  versus 
nourishment  sand  volume  added  for  three  different  grain  sizes.  The  common 
characteristics  are: 

A,  =  0.10  m‘'^ 

=  Sp  =  Jeq  =  1:10 

AA„,i,  =  0.2  m 
B  =  1.5  m 
A.  =  6  m 


These  results  are  in  the  form  of  Figure  6,  which  was  based  on  perfectly 
sorted  sediment.  For  the  upper  curve  in  Figure  17,  which  applies  for  sediment 
coarser  than  the  native,  there  is  a  transition  from  intersecting  to  non-intersect¬ 
ing  profiles  at  an  added  volume  of  600  mVm.  For  the  other  two  sediments, 
all  equilibrium  profiles  are  non-intersecting. 

Summary  and  Conclusions  for  Methodology 
Presented 


Based  on  earlier  work,  three  distinct  types  of  equilibrium  profiles  can  exist: 
non-intersecting,  intersecting,  and  submerged. 
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profile  and  nourishment  material  grain  size  distribution 


INTERSECTING  NON-INTERSECTING 
PROFILES  PROFILES 


Methods  have  been  developed  and  illustrated  with  examples  to  calculate  the 
non-intersecting  and  intersecting  equilibrium  beach  profiles  resulting  from 
beach  nourishment.  With  sand  of  different  characteristics,  the  method  can 
accommodate  varying  ranges  of  realism,  from  idealized  initial  profile  and 
grain  size  distribution  nourished  profiles  represented  by  analytical  forms,  to 
the  most  realistic  case,  in  which  the  initial  profile  and  nourishment  grain  size 
distribution  are  arbitrary  and  are  user-specified.  The  equilibrium-nourished 
profile  is  based  locally  on  the  differential  equation  for  equilibrium  for  a 
distribution  of  sizes.  Examples  are  presented  illustrating  the  influence  of 
various  parameters. 
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4  Laboratory  Studies 


Introduction  and  Description  of  Facilities 

An  exploratory  series  of  laboratory  studies  were  conducted  to  investigate 
profile  response  to  waves  of  different  characteristics  for  the  case  of  a  poorly 
sorted  sediment.  Six  experiments  were  conducted  using  the  same  sand,  but 
with  different  initial  beach  slopes  and  wave  conditions. 

The  test  conditions  for  the  laboratory  experiments  are  presented  in  Table  1, 
and  a  schematic  of  the  facilities  is  shown  in  Figure  18. 


Table  1 

Experintental  Conditiorta,  Wave  Tank  Tasts 

Run  No. 

Wave  Period 
T(sec) 

Wave  Height 
H(ctnl 

Water  Depth 
h(cm) 

Initial  Slope 

1 

1.25 

9.0 

22.5 

1:10.85 

2 

1.25 

9.0 

22.0 

1:5.74 

3 

1.25 

11.0 

22.5 

1:9.93 

4 

1.25 

11.0 

21.0 

1:13.94 

5 

1.25 

8.5 

20.0 

1:24.27 

6 

1.25 

9.0 

21.0 

1:14.33 

Objectives 

The  general  objectives  of  the  laboratory  studies  were; 

a.  To  document  the  evolution  and  sorting  with  time  of  initially  planar 
beach  profiles  and  poorly  sorted  sediments. 

b.  To  compare  experimental  profiles  and  sediment  size  distributions  with 
predict^  values  based  on  techniques  developed  in  this  study. 
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1 8.  Schematic  of  laboratory  facilities 


Experiment  Procedures 


The  six  tests  commenced  with  the  initially  planar  beach  slopes  presented 
in  Table  1 .  Prior  to  establishing  the  initial  profile,  the  sand  from  across  the 
profile  was  mixed  to  approximate  uniformity.  After  establishing  the  uniform 
slope,  sand  samples  at  four  or  five  different  locations  across  the  profile  were 
collected  for  later  analysis. 

The  desired  wave  conditions  were  established  and  the  profile  documented 
and  sand  sampling  repeated  at  1,  S,  10,  and  24  hr  (Experiments  1-3)  and  1,  6, 
12,  and  24  hr  (Experiments  4-6).  The  wave  heights  were  measured  visually 
and  the  location  and  height  of  the  breaking  waves  were  documented  several 
times  during  the  test. 

Test  Results 

The  results  of  the  test  program  are  described  below  for  each  of  the 
experiments. 

Experiment  1 

The  profile  evolution  from  a  planar  slope  of  1:10.85  is  presented  in  Fig¬ 
ure  19  for  the  initial  profile  and  profiles  at  1,5,  10,  and  24  hr  of  testing.  It 
is  seen  that  only  minor  changes  occurred  between  10  and  24  hr,  indicating 
that  the  system  had  approached  equilibrium.  The  general  characteristics  of  the 
final  profile  relative  to  the  initial  include  a  concave  upward  profile  with  most 
of  the  sand  transported  seaward  and  only  a  minor  amount  transported  land¬ 
ward  to  form  a  berm  feature.  Figure  20  presents  the  variation  with  time  of 
mean  grain  diameter  for  four  locations  across  the  beach  for  all  six  experiments 
at  0,  10,  and  24  hr.  For  Experiment  1  (Figure  20a),  it  is  clear  that  the  initial 
mean  grain  size  was  reasonably  uniform  across  the  beach  and  that,  with 
progressing  time,  the  coarser  sediments  were  transported  shoreward  and  the 
finer  sediments  seaward.  Figure  21  presents  the  initial  and  final  profiles  and 
the  grain  size  distributions  at  each  of  three  locations  across  the  profile,  and  the 
initial  grain  size  distribution.  Figure  22  shows  the  grain  size  distribution  at 
24  hr  at  six  locations  across  the  profile.  For  this  experiment,  substantial 
cross-profile  sorting  is  evident,  with  the  coarser  sediment  concentrated  near 
the  water  line. 

Experiment  2 

The  evolution  of  this  profile  is  presented  in  Figure  23  for  times  of  0,  1, 

5,  10,  and  24  hr.  This  profile  commenced  with  a  relatively  steep  uniform 
slope  (1:5.74)  and  practically  all  of  the  sediment  transport  was  seaward. 
Variations  with  time  of  mean  diameters  at  four  locations  across  the  profile  are 
in  Figure  20b.  There  is  substantially  less  pattern  to  the  mean  size 
distributions  compared  to  Experiment  1 .  There  has  been  some  reduction 
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Experiment  1 ,  measured  profiles  at  various  times 


ORICINni 


DISTANCE  (m) 

Experiment  1 ,  initial  and  final  profiles  and  initial  and  final  grain  size  distributions 
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'6  22.  Experiment  1 ,  grain  size  distributions  at  six  locations  across  the  profile  after  24  hr  of  testing 


ORIGINAL 


Experiment  2,  measured  profiles  at  various  times 


in  mean  grain  size  in  the  seaward  portion  of  the  profile,  and  some  size 
increases  toward  shore,  with  very  little  change  at  the  most  shoreward  location 
sampled.  The  initial  and  final  profiles  and  the  final  grain  size  distributions  at 
three  locations  across  the  profile  are  presented  in  Figure  24,  along  with  the 
initial  grain  size  distribution.  Figure  25  presents  the  grain  size  distribution  at 
six  locations  across  the  profile  after  24  hr  of  testing. 

Experiment  3 

Figure  26  presents  the  profile  evolution  for  Experiment  3  at  0,  1,5,  10, 
and  24  hr.  This  result  is  of  interest  in  that  wave  and  initial  slope  conditions 
are  nearly  the  same  as  presented  in  Experiment  1 ,  except  that  most  of  the 
sediment  transport  was  shoreward  in  Experiment  3.  Variations  with  time  of 
mean  grain  sizes  at  four  locations  across  the  beach  at  0,  5,  10,  and  24  hr  are 
presented  in  Figure  20c.  Figure  27  presents  the  initial  and  final  profiles  and 
the  grain  size  distributions  at  three  locations  across  the  profile.  Grain  size 
distributions  at  six  locations  across  the  profile  at  the  final  (24  hr)  survey  are 
presented  in  Figure  28.  In  general,  sorting  across  the  profile  has  occurred. 

Experiment  4 

Some  improvements  to  the  experiment  arrangement  were  made  before 
running  this  experiment,  including  water  level  control  and  recording  of  the 
waves.  The  conditions  involved  a  slope  of  1:13.94,  which  is  milder  than  in 
the  previous  experiments.  The  other  variables  were  kept  in  the  same  range. 

As  in  the  former  experiment,  irregular  waves  occurred. 

Inspection  of  Figure  29  shows  that  a  berm  and  a  bar  were  formed  as  in 
the  first  experiment.  The  volumes  do  not  match  as  well  as  in  previous  cases, 
which  is  due  to  lack  of  lateral  symmetry  and  the  substantial  consolidation 
observed  at  the  very  beginning  of  the  experiment.  Nevertheless,  erosion  at 
the  offshore  end  of  the  beach  was  observed.  Figure  29  shows  the  profiles  for 
U,  1,6,  12,  and  2‘1  hr  Grain  size  variations  are  shown  in  Figure  30.  The 
variations  with  time  of  mean  grain  sizes  at  four  locations  across  the  profile  are 
presented  in  Figure  20d.  Figure  31  provides  the  grain  size  distributions  at  six 
locations  across  the  profile  at  the  final  (24-hr)  survey. 

Although  a  clear  trend  in  grain  sizes  is  not  evident,  it  can  be  seen  that 
the  most  seaward  sample  has  about  the  same  distribution  as  the  original  and 
that,  this  time,  a  sorting  to  finer  sizes  has  been  achieved  at  the  berm. 

Experiment  5 

This  experiment  included  the  mildest  initial  slope  of  all  six  tests.  A  berm 
trapping  a  lagoon  was  formed  with  the  sand  apparently  originating  from  both 
onshore  and  offshore  sides  of  the  beach.  Also,  an  offshore  bar  was  formed 
clearly  with  sand  provided  by  the  zone  in  between.  The  evolution  of  the 
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Experiment  2,  initial  and  final  profiles  and  initial  and  final  grain  size  distributions 


Figure  25.  Experiment  2,  grain  size  distributions  at  six  locations  across  the  profile  after  24  hr  of  testing 


Experiment  3,  measured  profiles  at  various  times 
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Figure  27.  Experiment  3,  initial  and  final  profiles  and  initial  and  final  grain  size  distribution 
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Figure  28.  Experiment  3,  grain  size  distributions  at  six  locations  across  the  profile  after  24  hr  of  testing 
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Figure  29.  Experiment  4,  measured  profiles  at  various  times 
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Experiment  4,  initial  and  final  profiles  and  initial  and  final  grain  size  distributions 
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e  31 .  Experiment  4,  grain  size  distributions  at  six  locations  across  the  profile  after  24  hr 


beach  is  shown  in  Figure  32  for  0,  1,6,  12,  and  24  hr.  The  initially  mild 
slope  became  even  milder  on  the  average  below  the  still-water  level.  From 
Figure  33,  which  shows  the  grain  size  distributions  for  the  initial  and  final,  it 
can  be  concluded  that  the  concentration  of  fine  grain  sizes  offshore  is  greater. 
The  variation  with  time  of  mean  grain  sizes  at  four  locations  across  the  profile 
is  presented  as  Figure  20e.  Figure  34  provides  the  grain  size  distributions  at 
six  locations  across  the  profile  at  the  final  (24-hr)  survey. 

Experiment  6 

This  was  the  most  sophisticated  experiment  because  of  the  use  of  sand 
tracers.  The  evolution  of  the  beach  is  shown  in  Figure  35.  As  in  Experi¬ 
ment  1 ,  a  berm  was  formed  but  most  of  the  sand  was  transported  offshore 
where  it  formed  a  bar  feature. 

The  sand  size  distributions  for  the  initial  and  final  profiles  are  shown  in 
Figure  36.  It  is  clear  that  far  offshore  the  sand  is  finer  even  though  from  this 
figure  the  location  of  the  coarsest  sand  is  not  so  evident.  Some  shift  to 
coarser  sizes  from  the  initial  occurs  at  the  berm  and  at  3.0  m. 

The  tracers  were  followed  visually  as  carefully  as  possible,  but  this  was 
not  always  an  easy  task. 

The  initial  tracer  sieve  sizes  and  distributions  were  as  follows: 

a.  Blue:  #100  at  4.5  m. 

b.  Orange:  #100  at  shoreline. 

c.  Magenta:  #70  at  4.5  m. 

d.  Green;  #70  at  shoreline. 

e.  Yellow:  #50  at  4.5  m. 

/.  Red:  #50  at  shoreline. 

The  results  of  the  tracer  investigation  are  summarized  in  Table  2.  The 
following  can  be  concluded  from  the  analysis: 

a.  The  red  tracers  were  the  most  readily  tracked. 

b.  The  blue  tracer  was  lost  completely.  This  was  probably  due  to  the 
small  size  of  these  grains  (#100)  and  their  initial  location  in  relatively 
deep  water.  The  interpretation  is  that  when  the  action  of  the  waves 
began  they  were  suspended  and  spread  over  a  large  area  with  a  very 
low  concentration,  making  it  difficult  to  follow  their  path  or  even  to 
find  them  in  a  later  examination  of  the  sand. 
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Experiment  5,  measured  profiles  at  various  times 
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Experiment  5,  initial  and  final  profiles  and  initial  and  final  grain  size  distributions 
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Figure  34.  Experiment  5,  grain  size  distributions  at  six  locations  across 


Experiment  6,  measured  profiles  at  various  times 
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Figure  36.  Experiment  6,  initial  and  final  profiles  and  initial  and  final  grain  size  distributions 


Irable  2 

1  Results  of  the  Tracer  Analysis 

|Tim« 

Type  of 
Obaarvation 

Oranga 

Groan 

Rad 

YaVow 

Magenta 

1  HR 

Visual 

At  1.5  m 

Dyed  zona 

At  1.6-1. 7  m 

Dyed  zone 

At  1. 7-2.0  m 

Spots 

Sample  tubes 

At  1 .5  m 

On  the  top 

At  1.7  m 

Few  grams 

■ 

Accumulation 

tuba 

■ 

At  1 .9  m 

Many 

Many 

At  2.5  m 

Some 

Some 

Some 

At  3.0  n. 

At  3.5  m 

Some 

At  4.0  m 

Slightly  dyed 

miHi 

Many 

At  5.0  m 

Few 

Many 

1  1 

6  HR 

Visual 

Layers  by 
the  glass 

■ 

■ 

Accumulation 

tube 

■ 

At  2.0  m 

Few 

At  2.5  m 

Few 

At  2.9  m 

Some 

At  3.1  m 

Few 

Few 

At  3.2  m 

Few 

\  1 

1  24  HR 

Sample  tubes 

■ 

At  1.95  m 

Layer  at  4  cm 
below  top 

Few 

■ 

■ 

At  2.30  m 

Few  at  1-5  cm 
from  top 

■ 

1  (Continued) 
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Table  2  (Concluded)  | 

Tima 

Type  of 
Observation 

Orange 

Green 

Rad 

YoUow 

At  2.50  m 

Layer  at  1 .5 

cm 

below  top 

Few 

At  2.80  m 

Few  top 

Few  top 

At  3.00  m 

3  Grains 

At  3.50  m 

Few 

Few 

At  3.80  m 

At  4.20  m 

Few 

At  4.60  m 

Few 

At  5.10  m 

Layer  at 

3.5  cm  below 
top 

c.  Problems  similar  to  those  with  the  blue  tracer  could  have  occurred 
with  the  orange;  however,  as  the  placement  of  the  orange  was  closer 
to  the  shoreline,  the  area  over  which  it  was  spread  was  smaller,  and 
the  concentration  higher.  This  allowed  these  particles  to  be  identified 
upon  the  conclusion  of  the  experiment. 

d.  The  tracers  that  were  initially  located  at  the  shoreline  were  transported 
and  sorted  on  the  beach  face  during  the  first  hour  of  wave  action. 

They  remained  there  until  the  end  of  the  experiment.  At  that  time, 
they  were  located  4  cm  below  the  surface  corresponding  to  the  profile 
during  the  first  hour. 

e.  The  orange  tracer  was  found  at  the  top  of  the  berm  while  the  red,  not 
so  evenly  spread,  was  found  at  the  low  part  of  the  beach  face.  Recall 
that  the  orange  tracer  was  the  finest  and  that  the  red  was  the  coarsest. 

/.  At  the  bottom  of  the  profile,  just  before  the  offshore  bar,  tracers  of 
different  colors  were  found,  indicating  that  transport  has  occurred  in 
both  directions-seaward  and  landward-of  this  point. 

g.  No  other  evidence  was  found  to  indicate  major  patterns  in  sediment 
transport. 

Other  results  concerning  all  experiments  are  presented  below  and  will  be 
used  later  to  develop  additional  conclusions. 

Figures  20f  and  37  present  the  remaining  grain  size  information  for 
Experiment  6. 
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By  comparison,  the  profiles  equilibrated  after  12  hr  of  wave  action  and 
then,  only  approximately.  However,  the  main  features  of  the  final  profile 
were  established  very  fast;  that  is,  the  beach  face  slope,  the  sandbar,  and  the 
berm  were  formed  in  the  first  6  hr  and  after  ,hat,  they  only  migrated  or 
modified  their  volumes  slightly.  Layers  of  sand  were  added  or  removed  from 
the  existing  features,  but  the  slopes  that  dominate  the  main  characteristics  of 
the  beach  remained  almost  constant.  The  tracers  (Experiment  6)  onl^ 
confirmed  the  way  in  which  the  main  beach  features  developed.  At  the  very 
beginning  of  the  experiment,  the  waves  carried  the  finest  tracers  up  the  beach 
face.  Furthermore,  the  coarsest  tracer  type  was  spread  over  a  larger  area. 
After  some  hours,  a  berm  had  been  formed  and  the  beach  face  became 
steeper.  The  sediments  were  no  longer  able  to  be  transported  up  the  beach  as 
far  as  before,  so  the  orange  tracers,  which  were  transported  to  the  farthest 
onshore  location,  remained  uncovered  until  the  end  of  the  experiment.  On  the 
other  hand,  the  other  tracers  were  covered  by  successive  layers  being  found 
after  24  hr  at  some  depths  below  the  surface. 


Comparison  of  Laboratory  Data  with  Predictions 

In  this  section,  the  six  "final"  (i.e.  24-hr)  profiles  obtained  in  the  labora¬ 
tory  experiments  were  compared  with  computed  profiles  based  on  the  cross- 
shore  varying  mean  grain  sizes  as  documented  in  the  laboratory.  It  is  stressed 
that  the  comparisons  to  be  presented  are  "blindfolded"  in  the  sense  that  they 
do  not  incorporate  any  calibrations  or  adjustments  to  improve  the  fit.  Rather, 
the  computed  profiles  are  based  on  the  empirical  relationship  between  the 
sediment  scale  parameter  A,  and  the  sediment  size  D,  presented  in  Figure  2. 
Comparisons  are  presented  for:  (a)  profiles  for  parameterized  fits  to  the 
actual  A  versus  y  distributions,  and  (b)  profiles  for  the  local  A  values. 


Comparison  Based  on  Parameterized  Fit  to  the 
Actual  A  Values 

The  method  will  first  be  described,  then  the  results  of  the  comparison 
presented. 

The  procedure  considers  cross-shore  variations  of  A  of  the  forms 

A(y)  =  Ao  (9) 


A(y)  =  A„  *  my 


(10) 
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Figure  37.  Experiment  6,  grain  size  distributions  at  six  locations  across  the  profile  after  24  hr 


in  which  y  is  the  distance  offshore  from  the  still-water  line,  ^4.  is  the 
(idealized)  sediment  scale  parameter  evaluated  at  the  shoreline,  and  k  and  m 
are  empirical  coefficients  describing  the  best-fit  variation  in  the  cross-shore 
direction. 

The  general  expression  for  equilibrium  beach  profiles  is 

(11) 

dy 

in  which  A  is  the  local  value  of  the  sediment  scale  parameter.  Substituting 
Equation  9  into  Equation  1 1  and  integrating 

and  considering  the  variation  given  by  Equation  10, 

h(y)  =  [(Ao  ^  my)^^-Ai^ 

The  method  described  above  was  applied  as  follows.  The  distribution  of 
the  A  parameter  across  the  profile  was  determined  for  each  of  the  final 
sediment  samples  based  on  the  mean  diameter  and  transforming  to  A  via 
Figure  2.  Next,  the  best  least-squares  representations  of  these  data  by  Equa¬ 
tions  9  and  10  were  established.  Finally,  the  calculated  profiles  were  based 
on  Equations  12  and  13.  Figures  38-43  present  a  comparison  of  the  calculated 
and  measured  profiles  for  the  six  experiments.  These  results  are  discussed 
briefly  below. 

For  Experiments  1,5,  and  6,  the  predicted  profile  is  somewhat  steeper 
than  the  measured,  whereas  for  Experiments  2  and  3,  the  measured  is  steqier 
than  the  predicted  and  for  Experiment  4,  there  is  reasonable  agreement. 


Comparisons  Based  on  Localized  A  Values 

In  this  comparison,  which  was  only  applied  to  Experiments  1,  2,  and  3, 
the  calculated  profile  was  based  on  the  local  A(y)  value  as  given  by 

h(y+dy)  =  [h^(y)  +  A^(y)dyf^ 


which  follows  directly  from  the  differential  form  (Equation  11).  Results  of 
applying  Equation  14  to  calculate  profiles  are  presented  in  Figures  44  and  45. 
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Results  are  somewhat  similar  to  those  presented  and  described  for  the 
method  using  parameterized  fits  to  the  actual  A  versus  y  distributions  and  will 
not  be  discussed  further. 
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Experiment  1,  equilibrium  beach  profiles.  Top:  exponential  variation  of  A;  middle:  linear  variation  of  A;  bottom: 
Average  A 


CfiLCULftlEO  profile 
NEflSUREO  profile 


Experiment  2,  equilibrium  beach  profiles.  Top:  exponential  variation  of  A;  middle;  linear  variation  of  A;  bottom: 
Average  A 
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Experiment  3,  equilibrium  beach  profiles.  Top:  exponential  variation  of  A;  middle:  linear  variation  of  A;  bottom: 
Average  A 
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Experiment  5,  equilibrium  beach  profiles.  Top:  exponential  variation  of  A;  middle:  linear  variation  of  A;  bottom: 
Average  A 
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DISTANCE  (tn) 

Experiments  1 ,  2,  and  3,  comparison  of  predicted  and  measured  profiles  for  exponential  fit  to  A  values  as 
determined  from  measured  mean  sediment  sizes 
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5  Field  Data 


Two  sets  of  field  data  were  selected  for  comparison  and  evaluation  using 
the  general  methodology  presented.  Field  data  sites  were  Delray  Beach, 
Florida,  and  Jupiter  Island,  Florida;  both  sites  have  experienced  multiple 
renourishment  events.  As  with  most  field  results,  the  data  sets  are  not  as 
complete  as  desired.  However,  even  though  incomplete  and  complex,  field 
data  are  valuable,  as  they  contain  no  scale  effects.  The  following  sections 
present  field  data  from  these  two  sites  and  utilize  the  available  data  to  evaluate 
the  general  methodology  presented  here. 

Delray  Beach,  Florida 

As  shown  in  Figure  46,  Delray  Beach  was  nourished  in  1973,  1978,  and 
1984.  Figure  47  presents  the  northerly  and  southerly  iiinits  of  the  1984 
nourishment.  The  mean  diameter  of  this  nourishment  material  is  reported  to 
be  0. 16  mm,  whereas  that  of  the  native  material  is  0.22  mm. 

Figure  48  portrays,  for  Station  180.88,  the  1988  distribution  of  mean 
sediment  size  across  the  beach  profile  and  Figure  49  presents  the  associated 
A  parameter  with  an  exponential  fit.  The  same  information  is  presented  in 
Figures  50  and  51  for  Station  184.88  and  in  Figures  52  and  53  for  Station 
187.88.  Additional  profile  and  sediment  size  information  for  Delray  Beach  is 
contained  in  Appendix  C.  Two  types  of  profile  comparisons  were  carried  out 
as  described  below. 

The  first  type  of  profile  comparison  is  "blindfolded"  in  the  sense  that  only 
the  best-fit  exponential  A  relationships  (Figures  49,  51,  and  53)  were  used 
with  Equation  12  for  the  calculated  profiles.  These  comparisons  are  presented 
in  Figure  54.  In  general,  the  comparisons  are  considered  to  be  quite  good, 
with  the  deviations  for  depths  greater  than  4  to  5  m  believed  to  be  due  to 
nourishment  material  not  equilibrating  to  greater  depths. 

The  second  type  of  profile  comparison  is  also  based  on  exponential  distri¬ 
butions  of  A  values  using  best  least  squares  fit  to  the  measured  profile  data. 
These  results  are  presented  in  Figure  55  and,  again,  the  agreement  is 
generally  good. 
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Figure  46.  Nourishment  events  at  Delray  Beach,  Florida,  and  subsequent  volume  changes  (after  Coastal  Planning  and 
Engineering,  Inc.  (1988)1 


Figure  47.  Location  map  ot  Delray  Beach,  Florida,  nourishment  project  (after 
Coastal  Planning  and  Engineering,  Inc.  (1988)) 
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Average  grain  size  variation  across  Profile  180.88,  Delray  Beach,  Florida,  1988 
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Figure  52.  Average  grain  size  variation  across  Profile  187.88,  Delray  Beach,  Florida,  1988 
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Comparison  of  computed  and  measured  profiles.  Delray  Beach,  Florida,  1 988.  Computed  profiles  are  best  fit  based 
on  exponential  A  parameter 


Jupiter  Island,  Florida 


This  site  has  been  nourished  five  times,  with  the  first  nourishment  also 
occurring  in  1973.  The  nourishment  history  is  summarized  in  Table  3. 
Because  no  information  is  available  describing  the  grain  size  distribution 
across  the  beach,  comparisons  are  limited  to  information  available  in  profiles, 
both  pre-  and  post-nourishment.  The  nourished  area  and  the  profile  designa¬ 
tion  are  shown  in  Figure  56.  The  mean  diameter  of  the  nourishment  material 
is  reported  to  be  0. 12  nun  compared  to  the  native  of  0.20  mm. 


Table  3  I 

Summary  of  Nouriahmant  Hiatory  (Peat  19731  at  Juphar  bland  | 

Yaar 

Sagmant 

Sagmant 
Langth,  m 

Volume,  Cubic  Matara  g 

1973 

1 

8,120 

1,850,000 

1974 

1 

2,800 

750,000 

1  1 

1977 

1 

790 

200,000 

1977 

2 

1,100 

163,000 

1978 

1 

2,330 

6'>0.000 

1  1 

1983 

1 

1,780 

454,000 

1983 

2 

960 

311.000 

1  1 

1987 

1 

950 

287.000 

1987 

2 

3,300 

1.067,000 

1987 

3 

1,080 

353.000 

i  1 

Total 

6,085,000 

Figure  57  presents  best  least  squares  fits  to  the  12  available  1987  measured 
profiles  in  which  an  exponential  variation  of  A  with  offshore  distance  has  been 
utilized.  The  general  fit  to  the  measured  profiles  is  considered  good,  although 
at  some  profiles  there  is  a  nearshore  rock  reef  (e.g..  Profile  J  13.87, 

Figure  57,  sheet  2)  which  protrudes  above  the  sand  surface. 
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Figure  56.  Project  limits  and  profile  designations,  Jupiter  Island  project  (after 
Strock,  Arthur  V.,  and  Associates,  Inc.  1981)) 
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on  exponentiaf  /I  parameter  distribution  (Sheet  1  of  4) 
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distance 


The  second  type  of  analysis  presented  compares  the  effective  A  parameter 
in  the  shallower  and  deeper  portions  of  the  profiles  before  nourishment  (1973) 
and  after  nourishment  (1981  and  1987).  Based  on  the  exponential  distribution 
that  resulted  in  the  profile  fits  in  Figure  57,  the  A  values  near  the  shore  (A„), 
at  the  end  of  the  profile  lines,  and  the  average  of  these  two  are  compared  for 
1983,  1981,  and  1987  in  Figures  58,  59,  and  60,  respectively.  In  those  plots, 
north  is  to  the  right.  It  is  noted  that  for  pre-nourishment  conditions  (1973), 
there  was  a  substantial  difference  between  the  shallow-  water  and  end-of-Iine 
A  values.  In  1981  (Figure  59),  after  nourishment,  the  range  had  decreased 
some;  however,  the  mean  was  about  the  same.  This  is  interpreted  as  due  to 
the  beach  berm  being  displaced  by  the  relatively  finer  nourishment  sand  into 
deeper  water,  such  that  the  nearshore  and  end-of-line  A  values  were  more 
nearly  the  same.  Finally,  in  1987  (Figure  60)  the  A  values  at  the  shoreline 
and  end-of-line  are  nearly  the  same.  This  is  interpreted  as  the  nourished 
sediment  being  transported  over  the  entire  profile  such  that  the  A  parameter  is 
approximately  uniform  across  the  profile. 

A  simpler  type  of  analysis,  but  similar  to  that  presented  in  the  preceding 
paragraph,  was  carried  out  based  on  the  overall  characteristic  of  the  profile. 

In  particular,  if  h  =  Ay^'^  is  appropriate,  then  the  average  slope  sto  the  end  of 
line  where  the  water  depth  is  h'  at  an  offshore  distance  y'  is 

-  Af  (15) 


Hence  the  A  value  based  on  overall  slope  A,  is 

\  =  [5(/i')'T 


Similarly,  it  can  be  shown  readily  that  the  volume  per  unit  length  V  to  the 
end  of  the  line  is 


V  .  I 


3  (h')^'^ 
5  (Av)^ 


(17) 


Thus  the  A  value  based  on  volume  Ay  is 


Ay 


3(h')^ 
5  V 


TV3 


(18) 
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to  A  parameter  distribution.  Jupiter  Island,  Florida 
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In  interpreting  the  differences  between  A,  and  A^,  it  is  noted  that  if  the 
sediment  is  coarser  nearshore,  the  values  of  A  will  be  greater  than  offshore, 
and  it  is  to  be  expected  that  Ay  >  A„  as  will  be  seen.  Figures  61,  62,  and 
63  each  present  A  values  based  on  slopes  and  volumes,  and  the  average  of  the 
two  for  1973,  1981,  and  1987,  respectively.  The  results  and  interpretation 
are  generally  similar  to  those  presented  for  Figures  58,  59,  and  60.  Prior  to 
nourishment  (Figure  61),  the  A  values  based  on  volumes  were  somewhat 
greater  than  those  based  on  slopes  and  the  average  A  value  was  approximately 
0.10  Following  nourishment,  the  relative  difference  between  the  two  A 
values  is  approximately  the  same,  except  that  the  average  A  has  decreased 
somewhat.  Finally,  in  1987,  after  further  nourishment,  the  A  value  had 
decreased  further  to  less  than  0.08  m'^^  and  the  A  values  based  on  area  and 
volume  are  nearly  the  same. 

Conclusions  Based  on  Field  Data 

Based  on  analysis  of  available  field  data  from  two  nourishment  sites  at 
Delray  Beach,  Florida,  and  Jupiter  Island,  Florida,  the  following  are  con¬ 
cluded.  A  "blindfolded"  comparison  of  predicted  and  measured  (1988) 
profiles  for  Delray  Beach  shows  good  agreement  (Figure  54).  The  computed 
profiles  in  this  comparison  are  based  on  an  exponential  fit  to  the  A  values 
associated  with  the  mean  grain  sizes.  Differences  exist  primarily  for  depths 
greater  than  4  to  5  m  and  may  be  attributed  to  the  nourished  profile 
equilibrating  only  to  these  depths.  Additionally,  by  allowing  the  parameters 
in  the  exponential  A  representation  to  be  free,  a  goo'd  fit  is  obtained  across  the 
entire  profile  (Figure  55).  At  Jupiter  Island,  since  there  are  no  cross-shore 
sediment  size  data,  the  analysis  concentrated  on  the  variation  of  A  values  in 
the  nearshore  and  near  the  end  of  the  line.  Two  methods  were  employed  and 
although  the  general  results  of  the  two  methods  were  similar,  the  quantitative 
results  differed.  Based  on  an  exponential  A  fit,  prior  to  nourishment  the 
shoreline  A  value  was  substantially  greater  than  that  at  the  end  of  the  line.  By 
1987,  these  two  values  were  nearly  the  same  and  the  mean  value  had  not 
changed  appreciably  from  the  pre-nourished  values.  A  second  approach  to 
examining  the  A  values  nearshore  and  near  the  end  of  the  line  is  based  on  the 
average  slopes  and  volumes  of  the  measured  profiles.  Again,  prior  to 
nourishment,  the  two  A  values  differed,  whereas  in  1987,  they  were  virtually 
identical.  Results  using  the  second  method  differed  from  those  based  on  the 
exponential  A  fit  in  that  the  post-nourishment  average  A  value  decreased  to 
approximately  0.08  m"^  from  the  pre-nourishment  value  of  0.10  m’'^.  This 
result  is  qualitatively  consistent  with  the  use  of  nourishment  sediment  that  is 
finer  than  the  native  sediment. 
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Figure  62.  Longshore  variations  of  A  values  based  on  slopes,  volumes,  and  the  average  of  the  two,  Jupiter  Island,  Florida, 
1981 
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Figure  63.  Longshore  variations  of  A  values  based  on  slopes,  volumes,  and  the  average  of  the  two,  Jupiter  Island,  Florida, 
1987 


6  Summary,  Conclusions,  and 
Recommendations  for 
Further  Research 


Summary 


A  method  has  been  developed  and  illustrated  by  application  for  the 
prediction  of  the  two-dimensional  beach  profile  equilibrium  resulting  from  the 
placement  of  a  specified  volume  of  a  well-mixed  sand  of  arbitrary  size 
distribution.  It  is  shown  that  upon  placement  of  a  given  volume  of  material, 
three  types  of  profiles  can  result,  depending  on  the  material  size  characteris¬ 
tics,  the  volume,  berm  height,  and  closure  depth.  These  are:  (a)  intersecting, 
(b)  non-intersecting,  and  (c)  submerged  profiles.  The  method  is  applicable  to 
the  first  two  types.  The  method  assumes  that,  locally,  the  profile  is  in  equi¬ 
librium  with  the  profile  scale  characteristics  consistent  with  a  relationship 
developed  by  Moore  (1982).  An  iterative  method  is  employed  which  ensures 
that  the  volume  eroded  from  the  placed  profile  is  equal  to  that  deposited  sea¬ 
ward. 

Applications  provided  to  illustrate  the  method  have  included  idealized  grain 
size  distributions  and  profiles,  a  specified  grain  size  distribution  and  profile, 
and  a  range  of  various  mean  grain  sizes  and  sorting. 

Limited  small-scale  wave  tank  tests  were  conducted  to  investigate  sediment 
sorting  occurring  due  to  profile  evolution  from  an  initially  planar  slope. 
Surveys  and  surface  sand  samples  were  taken  across  the  profile  at  approxi¬ 
mately  0,  1,5,  10,  and  24  hr  after  conunencement  of  testing.  The  sand  sam¬ 
ples  were  later  analyzed  for  grain  size  distribution.  Although  the  results  were 
not  completely  consistent,  it  was  found  that  in  all  six  cases,  the  mean  sediment 
size  decreased  with  seaward  distance  from  the  equilibrium  shoreline. 

The  concepts  of  the  method  were  compared  where  possible  with  field  data 
from  Delray  Beach,  Florida,  and  Jupiter  Island,  Florida,  both  of  which  have 
been  nourished  on  multiple  occasions.  It  was  found  that  the  profile  shapes 
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could  be  predicted  reasonably  well  using  the  methods  employed  here.  Addi¬ 
tionally,  Ae  effects  of  nourishment  at  Jupiter  Island  caused  changes  in  profile 
shape  consistent  with  the  concepts  employed. 

Conclusions 

Previous  methods  for  assessing  the  relative  quality  of  sediment  sources  for 
beach  nourishment  have  compared  the  textural  properties  of  borrow  and  native 
material  and  result  in  an  "overfill  factor."  However,  these  methods  do  not 
provide  a  rational  capability  for  predicting  the  equilibrium  dry  beach  width, 
an  important  design  quantity  and  one  that  can  affect  significantly  the  project 
benefit/cost  ratio. 

Comparisons  of  predictions  with  laboratory-derived  profiles  were  reason¬ 
ably  encouraging,  with  the  predictions  overestimating  the  depths  in  three  ex¬ 
periments,  underestimating  the  depth  in  two  experiments,  and  generally  agree¬ 
ing  in  the  remaining  experiment. 

The  laboratory  cross-shore  mean  grain  sizes  seaward  of  the  equilibrium 
shoreline  were  generally  consistent  with  those  expected  and  incorporated  in 
the  method.  However,  the  data  included  a  considerable  distribution  of  sizes  at 
a  given  location,  whereas  the  method  employed  considered  the  particles  at  any 
location  to  be  perfectly  sorted. 

Blindfolded  tests  conducted  using  the  Delray  Beach,  Florida,  cross-shore 
grain  size  data  yielded  good  comparisons  between  measured  and  predicted 
beach  profiles.  The  agreement  was  poorer  for  water  depths  exceeding  4  to 
5  m,  beyond  which  it  is  believed  that  the  nourished  profile  either  did  not  ex¬ 
tend  or  did  not  equilibrate. 

The  method  presented  in  this  report  is  considered  to  be  one  step  toward  a 
rational  procedure  for  assessing  the  complex  performance  of  nourishment 
projects  with  realistic  sediment  characteristics. 

Recommendations  for  Further  Research 

A  great  deal  remains  to  be  accomplished  toward  the  objective  of  developing 
rational  procedures  for  predicting  the  performance  of  beach  nourishment 
projects  with  realistic  characteristics.  A future  beach  nourishment  projects 
should  be  monitored  carefully  and  completely.  Of  primary  importance  are: 

(a)  the  three-dimensional  evolution  with  time;  (b)  the  forcing  functions, 
i.e.,  the  directional  wave  spectrum,  and,  if  near  an  inlet,  the  currents;  and 
(c)  the  initial  and  evolving  grain-size  characteristics  in  both  the  cross-shore 
and  longshore  dimensions. 
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Methods  should  be  developed  for  predicting  the  cross-shore  distribution  of 
sediment  sizes.  There  is  a  need  to  develop  procedures  for  representing  the 
onshore  transport  and  sorting  of  sediment  placed  as  "profile  nourishment." 
i.e.,  as  an  underwater  deposit. 

Large-scale  wave  tank  tests  should  be  conducted  of  cross-shore  profile  and 
sediment  sorting  evolution.  Finally,  research  should  be  directed  toward  devel¬ 
opment  of  a  process-based  method  of  predicting  cross-shore  profile  and  sedi 
ment  sorting  evolution. 
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o  o  n  non 


C  THIS  PROGSAH  EQPR.FOR 

C  *************************'***************•*******************•********* 

C  *  THIS  PROGRAM  DEVELOPED  FOR  THE  COASTAL  ENGINEERING  RESEARCH  CENTER  * 
C  *  AND  CALCDLATES  THE  EQUILIBRIUM  BEACH  PROFILE  FOR  A  NOURISHMENT  * 
C  *  SAND  OF  SPECIFIED  SIZE  CHARACTERISTICS  October  23,1991  (Revised)  * 

C  *«***«*'****'**'*******«*'*«***********'***************'******'************** 

C 

COMMON/A/  y(800) ,HO(800) ,HP(800) ,HEQ(800) ,DC(800)  ,AC(800) , 

1  PCC(800) ,VUSEDC(800) ,P(40) ,D(40) , A(40) , DI (40) ,AI(40) 

DIMENSION  WORD(20) ,VGENV(50) ,VUSEDV(50) ,PSV(50)  ,DYEQV(50)  , 

1  1FLAG(6) 

OPEN (UNIT=5 , FILE= ' EQPR . INP' , STATUS= ' OLD ' ) 

OPEN (UN1T=6 , FILE= ' EQPR . OUT' , STATUS= ' NEW ' ) 

122  FORMAT (20A4) 

123  FORMAT(6F8.2,2I6) 

124  FORMAT(8F8.2) 

126  FORMAT(//,'  STARTING  LOOP  ITS  =  ',16,'  DYEQ  =  ',F8.3) 

160  FORMAT(8I6) 

162  FORMAT(F8.2,3I6,2F8.2) 

164  FORMAT(8I6) 

166  FORMAT(I6,3F8.2) 

167  F0RMAT(I6,4F10.3) 

168  FORMAT ('  IMAX  IFLAG2  lEND  HEND  DYEQO  DYEQ  VUSED' , 

1  5X,'  VGEN  HCROSS') 

170  FORMAT(3 (I6,2F8.2) ) 

178  FORMAT('  I  Y(I)  HO(I)  HP(I)  HEQ(I)  AC(1)',1X, 

1'  VDSEDC(I)  PCC(I)') 

180  FORMAT(I6,7F9.3) 

181  FORMAT('  ITS  =  ',13,'  DYEQ  =  ',F8.4,'  PSV  =  ',F8.3,'  VUSED  =  ' 
1,E10.4,  '  VGEN  =  ',E10.4,'  HO  =  ',F8.3,'  HP  =  ' , F8 . 3 , '  HEQ  =  ', 

2  F8.3) 

182  FORMAT ('  SOLUTION  REACHED,  NON-INTERSECTING  PROFILES') 

183  FORMAT('  SOLUTION  REACHED,  INTERSECTING  PROFILES') 

184  FORMAT (//) 

185  FORMAT('  I  DYEQ(I)  VGENV(I)  VUSED(I)', 

1,'  PSV(I)') 

C  NOTE:  ALL  DIMENSIONS  ARE  IN  METERS  * 

C  HO(I)=ORIGINAL  PROFILE  * 

C  HP{I)=PLACED  (NOURISHED)  PROFILE  * 

C  HEQ(I)=EQDILIBRIUM  PROFILE  * 

C  ************************************************************************** 

C  ESTABLISH  INITIAL  AND  PLACED  PROFILES 

C 

READ(5,122) (WORD(I) ,1=1,20) 

WRITE(6,122)  (WORDd)  ,1=1,20) 

WRITE (*,122) (WORD(I) ,1=1,20) 

CALL  PROCHAR ( JPTYPE , HO , HP , SP , VADD , IMAX , HBERM , DY , Y , DYO , HSTAR , 

1  SEQ,HMIX) 

READ  IN  CHARACTERISTICS  OF  PLACED  SAND 

CALL  SEDCHAR(LMAX,LIMAX,D,P,DI,AI) 

SUBROUTINE  ATRANS  TRANSFORMS  SEDIMENT  DIAMETERS  TO  A  VALUES 
CALL  ATRANS(LMAX,LIMAX,D,P,A,DI,AI) 

OUTER  HSTAR  LOOP 

READ (5,167) ITSMAX , DYEQ 1 
DDYEQ=DYEQl/20.0 
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IF ( DDYEQ . LT . 1 . 0 ) DDYEQ=1 . 0 

ITERMX=3 

PS=0.7 

VAVAIL-PS*VADD 

DYEQ-DYEQl 

C  THIS  LCWP  SOLVES  FOR  PS  FOR  TRIAL  DYEQ 

DO  200  ITS-^l.ITSHAX 
VGEN-VAVAIL 
VUSED=VAVAIL 
DYEQO=DYEQ 

DO  198  ITER=1,ITERMX 
DO  196  1-1,800 
AC(I)-0.0 
VUSEDC(I)=0.0 
PCC(I)=0.0 
196  HEQ(I)=0.0 
VAVAIL-VGEN 

CALL  PSSR(DYEQ,DY,HBERM,VAVAIL,VDSED, 

1  SEQ,IMAX,I21AX,HMIX,HSTAR,VGEN,ITS,ITSMAX, 

2  IM,IFLAG) 

C  WRITE ( * , * ) ITS , 1 PS , DYEQ , PS , VDSED , VGEN 

198  CONTINUE 
DYEQV( ITS) -DYEQ 
VGENV( ITS) -VGEN 
VUSEDV ( ITS ) -VUSED 
PSV ( ITS ) -VGEN/VADD 
PSCOR-PSV(ITS) 

VAVAIL=PSCUR*VADD 

WRITE (*,181) ITS, DYEQ, PSV (ITS) , VUSED, VGEN, HO (IM) ,HP(IM) ,HEQ(IM) 
IF(ITS.GT.l)GO  TO  199 
AA-VGEN-VUSED 
DDYEQ-DDYEQ*AA/ABS (AA) 

WRITE ( * , * ) ITS , DYEQ , DDYEQ , AA 
GO  TO  200 

199  CONTINUE 

AAC-VGENV ( ITS ) -VUSEDV ( ITS ) 

AAP=VGENV(ITS-1) -VUSEDV(ITS-l) 

IF(AAC.EQ.O.O)GO  TO  202 
IF (AAC/AAP.  LT.  0.0)  DDYEQ— DDYEQ/2 . 0 

200  DYEQ-DYEQ+DDYEQ 
202  CONTINUE 

WRITE(*,*) (I,IFLAG(I) ,1=1,6) 

4  00  CONTIN  j. 

800  CONTINUE 

WRITE (6, 184) 

WRITE (6, 178) 

WRITE (6, 180) (I,Y(I) ,HO(I) ,HP(I) ,HEQ(I) ,AC(I) ,VUSEDC(I) ,PCC(I) , 
1  I=1,IMAX) 

WRITE (6, 184) 

WRITE(6,185) 

WRITE(6,167) (I,DYEQV(I) ,VGENV(I) , VUSEDV (I) , PSV (I) ,I=1,ITSMAX) 
WRITE(6,184) 

IF(IFLAG(5) .EQ.1)WRITE(*,182) 

IF(IFLAG(6) .EQ.1)WRITE(*,183) 

IF(IFIAG(5) .EQ.1)WRITE(6,182) 

IF(IFLAG(6) .EQ.1)WRITE(6,183) 

CLOSE (UNIT-5) 

CLOSE (UNIT-6) 

STOP 

END 
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noon  n  n  n  n  n 


THIS  SUBROUTIHE  CALCULATES  THE  HEW  EQUILIBRIUM  PROFILE  AND 
DETERMINES  THE  TYPE  OF  PROFILE  (INTERSECTING  OR  NON-INTER) 


SUBROUTINE  PSSR(DYEQ, DY.HBERM, VAVAIL.VUSED, 

1  SEQ,IMAX,UfAX,HMIX,HSTAR,VGEN, 

2  ITS,ITSMAX,IM,IFLAG) 


COMMON/A/  Y(800) ,H0(800) ,HP(800} ,H£Q(800) ,DC(800) ,AC(800) , 

1  PCC(800) ,VUSEDC(800) ,P(40) ,D(40) , A(40) ,DI(40) ,AI(40) 

DIMENSION  HDP(5) ,HDC(5) ,IFLAG(6) 

20  FORMAT('  OPUT  SR  PSSR  , ITS-', 14,'  ITSMAX=',I4, 

1  '  DYEQ=' ,F8.3, '  VUSED-' , F8 . 3 , '  VGEN=',F8.3) 

22  FORMAT ( 315, 7F10. 3) 

23  FORMAT ( 315, 12F10. 3) 

24  FORMAT {'  REACHED  CROSS-OVER  POINT') 

160  FORMAT (8F8. 2) 

Il=DYEQ/DY+2 
DHP-0.0 

IFIAG(1)-1  MEANS  HP (I+l) >HEQ(I+1) 

IFLAG(1)=0 

IFLAG(2)=1  MEANS  INTERSECTING  PROFILES 
IFLAG(2)=0 

IFIAG(3)=1  MEANS  PROFILE  HAS  REACHED  CRITERION  FOR  END  OF 
INTERSECTING  OR  NON-INTERSECTING  PROFILE 
IFLAG(3)=0 

C  IFLAG(4)=1  MEANS  SEAWARD  OF  PLACED  PROFILE 

IFLAG(4)=0 

C  IFLAG(5)=1  MEANS  PROFILE  HAS  REACHED  CLOSURE  DEPTH 

IFLAG(5)=0 

C  IFLAG(6)=1  MEANS  INTERSECTION  OF  ORIGINAL  AND  EQUILIBRIUM  PROFILES 

IFLAG(6)=0 
HEQ(I1)=HBERM 
HMIXP-HMIX 
HCROSS=0 . 0 
IMM1=IMAX-1 
DYA=y(Il) -DYEQ 

VGEN=(DYA+HMIX/SEQ) * (DYA*SEQ+HMIX)/2 . 0 

VUSED=VGEN-DYA**2*SEQ/2 . 0 

WRITE ( * , * ) I 1 , DYA , SEQ , HMIX , VGEN , VUSED 

HDC(1)=HP(I1) 

HDC(2)=HEQ(I1) 

HDC ( 3 ) =HEQ ( I 1 ) +HMIX 
HDC(4)=HO(Il) 

YC=Y(I1) 

DO  100  1=11, IMMl 
IM=I+1 
DO  96  L=l,4 
96  HDP(L)=HDC(L) 

HMIXP=HMIXC 

Dyy=Dy 

YP=YC 

IEND=I+1 

HMIXC=HMIX 

YC=Y(I+1) 

C  FOLLOWING  SR  CALCULATES  NEW  EQUIL  PROFILE 

CALL  YNEW(I,Y,HEQ,VUSED,SEQ,VAVAIL,A,P,LMAX, 

1  DYEQ,HBERM,AC,VUSEDC,PCC) 

C  WRITE(6,*)I.HEQ(I+1) ,HO(I+l) 

HEND=HEQ(I+1) 
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C  IF(ITS.NE.8)GO  TO  36 

C  WRITE(6,23)I,IFIAG{1) ,1FIAG(4)  ,Y{I+1) ,HEQ(I+1) ,HO(I+l) ,HP(I+1) , 

C  1  VGEN,VOSED,VAVAIL,DHP,DYY,yp 

36  DH-HO(I+1)-HEQ(Ih1) 

IF(DH.LT.HMIX.AND.IFLAG(1)  .EQ-0)HMIXODH 
HDC(1)-HP(I+1) 

HDC(2)-HEQ(I+1) 

HOC ( 3 ) -HEQ ( 1+1 ) +HMIXC 
HDC(4)>«HO(I+l) 

IA=2 

IB=3 

IF(IFLAG(1) .EQ.O)GO  TO  40 
IB=1 

IF(IFLAG(4) .EQ.O)GO  TO  40 
IB=4 

IF(IFLAG(5) .EQ.1.0R.IFLAG(6) .EQ.l)GO  TO  200 
40  CALL  V0L(HDP,HDC,IA,1B,YP,YC,SEQ,HMIXP,VGEN,VUSED,IFLAG,HSTAR) 
100  CONTINUE 

WRITE  (*,20)  ITS ,  ITSMAX ,  DYEQ ,  VUSED ,  VGEN 
200  RETURN 
END 

C  - 

SUBROUTINE  PROCAL ( Y ,  I , AO ,  SO , HO, HBERM) 

C  - 

DIMENSION  Y(aOO) ,HO(800) 

20  FORMAT (216, 8F8. 3) 

IF(Y(I) .GE.-HBERM/SO)  GO  TO  40 
HO(I)=HBERM+(Y(I)-Y(l))*SO 
GO  TO  42 
40  H0(l)=0.2 
DO  80  K=l,20 

C  WRITE(*,20)I,K,Y(I)  ,HO(I)  , SO, AO, HBERM 

EPS=Y ( I ) +HBERM/SO-HO ( I ) /SO- ( HO  ( I ) /AO) **1.5 
DEPSDH=-1.0/SO-1.5*HO(I)**O.S/AO**1.5 
80  HO(I)=HO(I)-EPS/DEPSDH 
42  RETURN 
END 
C 

C  THIS  SUBROUTINE  CALLED  FOR  EACH  Y  VALUE  TO  DETERMINE 

C  HEQ(I+1) 

C  - 

SUBROUTINE  YNEW  ( I ,  Y ,  HEQ ,  VUSED ,  SEQ  ,  VAVAIL ,  A ,  P ,  IJIAX , 

1  DYEQ,  HBERM,  AC,  VOSEDC,PCC) 

C  - 

DIMENSION  Y(800) ,HEQ(800)  ,A(40)  ,P(40) ,AC(800) ,VUSEDC(800) , 

1  PCC(800) 

20  FORMAT(2I6,6F8.3) 

22  FORMAT ('  OUTPUT  FROM  SR  YNEW') 

C  NOTE  HBERM  IS  NEGATIVE 

C  WRITE(*,22) 

Y2ERO=-HBERM/SEQ+DYEQ 
IF(Y(I+1) .LE.YZERO)GO  TO  40 

C  THIS  SECTION  FOR  PORTION  OF  PROFILE  BELOW  WATER 

DYY=Y(I+1)-Y(I) 

HPREV=HEQ(I) 

DYY2=Y(I+l)-YZERO 
IF(DyY2.GE.DYY)GO  TO  34 

C  THIS  SECTION  FOR  FIRST  BELOW  WATER  POINT 

DYY=DYY2 
HPREV=0 . 0 
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non 


34  PC=VUSED/VAVAIL 
VUSEDC(I)=VUSED 
IF(PC.GT.1.0)PC=1.0 
PCC(I)=PC 

IF (PC. GT. 1.0) PC-1.0 

C  THIS  SUBROUTINE  DETERMINES  CURRENT  A  PARAMETER  VALUE 

PC=1.0-PC 

CALL  INTERP(P,A,PC,AB,IMAX) 

AC(I)=AB 

C  WRITE (6, 20)1, LMAX , VA VAI L, VUSED , VGEN , AB 

HEQ ( I+l) -HPREV+0 . 02 
DO  36  LL=1,10 
HBAR=0.5*(HPREV+HEQ(I+1) ) 
BB=1.0/SEQ+1.5*SQRT(HBAR)/AC(I)**1.5 
DHDy=1.0/BB 

HEQ ( I + 1 ) =HPREV+DHDY  *  DY Y 
36  CONTINUE 
GO  TO  102 

C  NEXT  STATEMENT  FOR  PORTION  OF  PROFILE  ABOVE  WATER 

40  HEQ(I+1)=HBERM+{Y(I+1)-DYEQ)*SEQ 
102  RETURN 
END 

C  - 

SUBROUTINE  INTERP(X, Y , XC, YC,NAS) 


DIMENSION  X(40),Y(40) 

22  FORMAT ('  OUTPUT  FROM  INTERP' ) 

24  FORMAT(2I6,6F8.3) 

26  FORMAT('  ******♦', 216 , 6F8 . 3) 

C  WRITE(*,22) 

NASM1=NAS-1 
DO  40  N=1,NASM1 
NC=N 

IF(XC.GE.X(N) .AND.XC.LE.X(N+1) .AND.X{N+1) .GT.X(N) )  GO  TO  42 
IF(XC.LE.X(N) .AND.XC.GE.X(N+1) .AND.X(N+1) .LT.X(N) )  GO  TO  42 
40  CONTINUE 
42  DX=X(NC+1) -X(NC) 

DY=Y(NC+1) -Y(NC) 

YC=Y (NC) + (XC-X (NC) ) *DY/DX 

WRITE{8,26)NC,NAS,XC,YC,DY,DX,X(NC) ,Y(NC) 

WRITE (*,22) 

WRITE(*, 24)NC,NAS,XC,YC 
RETURN 
END 
C 

C  THIS  SUBROUTINE  ESTABLISHES  THE  CHARACTERISTICS  OF  THE 

C  ORIGINAL  AND  PLACED  PROFILES 

C  - 

SUBROUTINE  PROCHAR ( JPTYPE , HO , HP , SP , VADD , IMAX , HBERM , 

1  DY,Y,DYO,HSTAR,SEQ,HMIX) 


DIMENSION  Y(800) , HO(800) , yV(lOO) , HVO (100) ,HP(800) 

20  FORMAT ( 216, 9F7. 2) 

22  FORMAT ('  ORIGINAL  PROFILE') 

24  FORMAT ('  PLACED  PROFILE') 

26  FORMAT(3(I6,2F8.2)  ) 

28  FORMAT('  ITERATION-' , 16 , 'IL=' 16 , '  VADDED-' , F8 . 2 , '  VUSED- ' , F8 . 2  , 
1  2X, 'DYO-' ,F8.3) 

32  FORMAT('  VADDED  =  ',F8.2,'  VUSED  =  ',F8.3,'  DYO  =  ',F8.3) 

31  FORMAT(6I6) 
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30  FORMAT (' CHECK' ,16, 8F8. 3) 

33  FORMAT (8F8. 2) 

READ  ( 5 , 2 0 )  JPTYPE ,  IMAX ,  DY ,  AO ,  SO ,  SP ,  VADD,  SEQ ,  HBERM ,  HSTAR ,  HMIX 
WRITE  ( * ,  2 0 )  JPTYPE ,  IMAX ,  DY ,  AO ,  SO ,  SP,  VADD ,  SEQ ,  HBERM ,  HSTAR ,  HMIX 
Y(1)=0.0 
DO  10  1=2, IMAX 
10  Y(I)=Y(I-1)+DY 

C  JTYPE=1,READ  IN  HO  (ORIGINAL  DEPTHS);  JTYPE=0 ,  CALCULATE 

C  DEPTHS 

IF ( JPTYPE. EQ.l)  GO  TO  42 
DO  40  1=1, IMAX 

40  CALL  PROCAL(Y,I,AO,SO,HO,HBERM) 

GO  TO  44 

42  READ(5,31)IMAXP 

READ(5,33) (YV(I) ,HVO(I) ,I=1,IMAXP) 

H0(1)=HV0{1) 

DO  43  1=2, IMAX 
XC=Y(I) 

CALL  INTERP(YV,HVO,XC, YC,IMAXP) 

43  HO(I)=YC 

44  CONTINUE 

C  WRITE(6,22) 

C  WRITE(6,26) (I,Y(I) ,HO(I) ,1=1, IMAX) 

GO  TO  60 

C  THIS  SECTION  ESTABLISHES  PLACED  PROFILE 

C  60  WRITE (6, 20)  IMAX, JPTYPE, VADD 
60  CONTINUE 
DYO=20.0 
DO  80  K=l,20 
DO  50  1=1, IMAX 
50  HP(I)=0.0 
VUSED=0.0 

C  WRITE (6, 20) IMAX, JPTYPE, VADD 

HP(1)=HBERM 
DO  78  1=2, IMAX 
IM=I 
AB=0 . 0 
IL=I 

IF(Y(I) .LE.DYO)HP(I)=HBERM 

IF(Y(I) .GT.DYO)HP(I)=HBERM+(Y(I)-DYO)*SP 

DHC=HO ( I ) -HP ( I ) 

DYY=DY 

AC=HO(I-l)-HP(I-l) 

IF(HP(I) .LT.HO(I) )  GO  TO  64 
AA=Y(I)-Y(I-1) 

AB=HO(I)-HP(I) 

AC=HO(I-l)-HP(I-l) 

DYY=  AC*DY/ (AC-AB) 

DHC=0 . 0 

64  VUSED=VUSED+0.5*DYY*(DHC+AC) 

C  WRITE(6,20)K,I,VADD,VUSED,AA,AB,AC,DYY 

IF(AB.LT.0.0)GO  TO  79 

78  CONTINUE 

79  DYO=DYO+ (VADD-VUSED) / (HO (IL) -HBERM) 

C  WRITE(6,30)IL,DYO,VADD,VUSED,HO(IL) , HBERM 
WRITE ( * , 2 8 ) K , I L , VADD , VUSED , DYO 

80  CONTINUE 

C  DO  82  I=IM,IMAX 
C  82  HP(I)=HO(I) 

C  WRITE(6,24) 
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C  WRITE(6,26) {I,Y(I) ,HP{I) ,1-1,IMAX) 
RETURN 
END 


SUBROUTINE  SEOC31AR(LKAX,I.ZMAX,0,P,DZ,AI} 


DIMENSION  0(40) , P(40) , 01  (40) ,AI(40} ,PD(40) 

20  FORMAT (16) 

21  FORMAT(2I6.6F8.2) 

22  FORMAT (8F8. 3) 

24  FORMAT('  OUTPUT  FROM  SUBROUTINE  SEDCHAR,  UlAX  «  ',16) 
26  FORMAT('  INPUT  PAIRS  OF  DI,AI  ') 

28  FORMAT('  INPUT  PAIRS  OF  D(liiai),P') 

29  FORMAT(I6,4F8.3) 

30  FORMAT (//) 

READ ( 5 , 2 1 ) USED , LMAX , XMU , SIG 
C  WRITE(6,21)JSED,LMAX,XMU,SIG 
LIMAX-19 
LMMl-IMAX-1 
AA-0.3989 

REAO(S,22)  (DI(L)  ,AI(L)  ,1^1, UMAX) 

WRITE(6,30) 

WRITE (6, 26) 

WR1TE(6,22)  (DI(L)  ,AI(L)  , 0=1, UMAX) 

IF(JSED.EQ.l)GO  TO  8 
DO  10  L-1,1MAX 
XL=L 

DD-0.6*SIG 

D(L) -XMD-3 . 0*SIG+ (XL-1) *DD 
BB- ( (D (L) -XMU) /SIG) **2/2 . 0 
PD ( L) -AA/SIG*EXP ( -BB) 

DC=*D(L) 

10  D(L)=2.0**(-DC) 

C  WRITE(6,22) (D(L) ,PD(L) ,L=1,11) 

P(l)=1.0 

P(11)=0.0 

DO  12  I/=2,LMM1 

12  P(L)=P(L-l)-(PD(L-l)+PD(L))*DD/2.0 
GO  TO  80 

8  READ(5,22) (D(L) ,P(L) ,L=1,IMAX) 

C  WRITE (6, 26) 

80  WRITE(6,30) 

WRITE(6,28) 

WRITE(6,22) (D(L) ,P(L) ,L=1,LMAX) 

RETURN 

END 

C 

C  - 

SUBROUTINE  ADET { I , A , P , AC , VUSED , VAVAI L ,  IMAX ) 


DIMENSION  A(40) ,P{40) ,AC(800) 

PC=VUSED/VAVAIL 

DO  20  L=2,LMAX 

LC=IJ1AX-L+1 

IF(PC.GT.P(LC) . OR.PC.lt. P(LC+1))  GO  TO  20 
DP=P(LC)-P{LC+1) 

A1=PC-P(LC+1) 

DA=A(LC)-A(LC+1) 

AC ( I ) =A ( LC+ 1 ) -DP*  DA/DP 
GO  TO  22 
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20  CONTINUE 
22  CONTINUE 
RETURN 
END 

THIS  SUBROUTINE  TRANSFERS  SEDIMENT  DIAMETER  ,D,  TO 
SEDIMENT  SCALE  FACTOR, A 


SUBROUTINE  ATRANS(LMAX,LIMAX,0,F,A,DI,AI) 


DIMENSION  D(40) ,P(40) ,A(40) ,DI(40) .AI(40) 
22  FORMAT  ('  OUTPUT  FROM  SR  ATRANS') 

24  FORMAT(8F8.3) 

26  FORMAT(I6,3F8.3) 

28  FORMAT (216, 5F8. 3) 

WRITE(6,28)IMAX 

WRITE(*,22) 

DO  20  L-1,LMAX 
DOD(L) 

WRITE  (  6 , 2  8 )  L ,  IMAX ,  DC 

CALL  INTERP(DI,AI,DC,AC,LIMAX) 

WRITE{6,26)L,DC,AC,P(L) 

20  A(L)*AC 

WRITE(*,22) 

WRITE(6,22) 

WRITE(6,24)  (A(L)  ,P(L) ,L=1,IMAX) 

RETURN 

END 


SUBROUTINE  VOL ( HOP , HDC , lA , I B , YP , YC , SEQ , HMIXP , 
1  VGEN,VUSED,IFLAG,HSTAR) 


DIMENSION  HDP(5) ,HDC(5) ,IFLAG(6) 

20  FORMAT('  IFLAG(4)  =  ',16) 

22  FORMAT(2I6,7F8,3) 

24  FORMAT('  IFLAG(l)  =  ',16) 

26  FORMAT('  IFLAG(5)  =',I6,'  IFLAG(6)  =',I6) 

HMIX=HMIXP 

IF(IA.EQ.2.AND.IB.EQ.3)GO  TO  62 
IF(IA.EQ.2.AND.IB.EQ.l)GO  TO  64 
IF(IA.EQ.2.AND.IB.EQ.4)GO  TO  66 
THIS  SECTION  FOR  GENERATION  OF  PRIMARY  VOLUME 
62  IF(HDP(1)  .LT.HDP(2) .AND.HDC(l) .GT. HDC(2 ) ) IFLAG (1) =1 
IF(IFLAG(1)  .EQ.O)GO  TO  300 

IF  ABOVE  CONDITION  IS  MET,  HAVE  REACHED  CROSS-OVER  POINT 
IA=1 
IB=2 

CALL  CROSS  ( HDP ,  HDC ,  lA ,  IB ,  YP ,  YC ,  YINT ,  HINT) 

DVMIX= (HDP ( 3 ) -HDP ( 2 ) ) * ( YINT-YP) +HMIX*  *  2/ ( 2 . 0*SEQ ) 

DVGEN= (HDP(2)-HDP(l))*(YINT-YP)/2.0 
VUSED=VDSED+DVMIX+  ( HDC  ( 1 )  -HDC  ( 2  )  )  *  (  YC-YINT)  /  2 . 0 
VGEN=VGEN+DVGEN+DVMIX 
IF(HDC(1)  -EQ.0.0)IFLAG(4)=1 
WRITE ( 6 , 2 4 )  I FLAG ( 1 ) 

WRITE ( 6 , 22 )  lA , IB , DVMIX , DVGEN , VUSED . VGEN , YINT , HMIX , SEQ 
GO  TO  400 

C  TO  HERE  IF  PAST  CROSS-OVER  POINT  BUT  BEFORE  PROFILE  INTERSECTION 

64  IF(HDP(1)  .LT.HDP(4) .AND.HDC(l) .GT.HDC(4) ) IFLAG(4)=1 
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c 

c 


c 

c 


c 


c 

c 

c 

c 


c 

c 

c 

c 


c 

c 


c 

c 

c 


c 


c 

c 

c 


IF(HDC(1) .EQ.0)IFIAG{4)=1 
IF(IFLAG(4) .EQ.O)GO  TO  300 
WRITE ( 6 , 20)  IFLAG (4 ) 

IA=4 

CALL  CROSS (HDP,HDC,IA. IB, yP,yC, VINT, HINT) 

VUSEI>=VDSED+ (HDP  ( 1 ) -HOP  ( 2  ) +HIX:  ( 1 ) -HDC  { 2  )  )  *  ( y C-y P)  / 2 . 0 
1  »(HDC(l)-HDC(4))*(yC-yiNT)/2-0 

WRITE(6,22)IA,IB,VOSED,yC,yP,yiHT,HDP(l) ,HDP(2)  ,HDC(1) 
WR1TE(6,22)IA,IB,HDC(2) ,HDC(4) 

GO  TO  400 

66  IF(HDP(4)  .GT.HDP{2)  .AND-HDC{4)  .LT.HDC(2)  )  IFLAG(6)=1 
IF(HDC(2) .GT.HSTAR)IFLAG(5)=1 

IF(IFLAG(5)  .NE.1.AND.IFLAG(6)  .NE.DGO  TO  300 
XM2=(HDC(2)-HDP(2) )/(yC-yP) 

XM4=(HDC(4)-HDP(4) )/ iVC-VP) 

IF(IFLAG(5)  .EQ.l)GO  TO  68 
yiNT=y P+ ( HOP ( 2 ) -HOP ( 4 ) ) / ( XM4 -XM2 ) 

GO  TO  70 

TO  HERE  IF  REACHED  EQUILIBRIUM  DEPTH 
68  yiNT=yP+{HSTAR-HDP(2) )/XM2 

VUSED=VUS£CH-  (HDP ( 4 ) -HDP ( 2 ) ) * (VINT-YP) 

1  +(XM4-XM2)*(YINT-YP) **2/2.0 

GO  TO  400 

7 0  DVUSED= (HDP ( 4 ) -HDP (21)* ( YINT-VP)  + 

1  (XM4-XM2)*(yiNT-yP)/2.0 

VUSEI>=VUSED+DVUSED 
WRITE(6,26)IFLAG(5) ,IFLAG(6) 

WR1TE(6,22)IFLAG(5) ,IFLAG{6) ,yP,yC,HDP(2) ,HDP(4)  ,HDC(2) , 

1  HDC(4) 

WRITE ( 6 , 2 2 )  lA , IB , XM2 , XM4 , TINT , DVUSED ,  VUSED 
GO  TO  400 

300  CALL  DVOL(HDP,HDC,IA,IB,yp,yC,DW) 

IF(ITS.EQ.8)WRITE(8,*)ITS,I,IA,IB,YP,yC,DW,HDP(2)  , 

1  HDP(3) ,HDC(2) ,HDC(3) 

IF(IFLAG(4)  .EQ.1)WRITE(6,22)IA,IB,VGEN,DW,HDP(2)  ,HDP(4)  , 
1  HDC(2) ,HDC(4) ,YP 
VUSED=VUSED+DVV 
IF(IFLAG(1) .GT.O)GO  TO  400 
IA=1 
IB=3 

CALL  DVOL(HDP,HDC,IA,IB,YP,YC,DW) 

IF(ITS.EQ.8)WRITE(8,*)ITS,I,IA,IB,YP,yC,DW,VGEN, 

1  HDP(l) ,HDP(3) ,HDC(1) ,HDC(3) 

VGEN=VGEN+DW 
400  RETURN 
END 


SUBROUTINE  DVOL(HDP, HDC,  lA,  IB,  YP,  YC,  DW) 


DIMENSION  HDC(5) ,HDP(5) 

20  FORMAT (216, 7F8. 2) 

DW=0 . 5*  ( YC-yP)  *  (HDP  ( IB)  -HDP (lA)  +HDC  (IB)  -HDC  ( lA)  ) 
WRITE(7,20)IA,IB,YP,YC,HDP(IA) ,HDP(IB) ,HDC(IA)  ,HDC(IB) ,DW 
RETURN 
END 


SUBROUTINE  CROSS (HDP , HDC ,IA,IB,YP,YC,yiNT,HINT) 
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DIMENSION  HDP(5) ,HDC(5) 

20  FORMAT (2 16, 8F8. 2) 

C  WRITE(6,20)IA,IB,YP,YC,HDP(IA) ,HDP(IB) ,HDC(IA)  ,HDC(IB) 

XMA=(HDC(IA)-HDP(1A) )/(YC-YP) 

XMB= (HDC ( IB) -HDP ( IB) ) / ( YC-YP) 

YINT=YP+(HDP(IA)-HDP(IB) )/(XMB-XMA) 

HINT-HDP(IA) +XMA* (YINT-YP) 

C  WRITE{6,20)IA,IB,XMA,XMB,YINT,HINT 

RETURN 
END 
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EXAMPLE  1,  CERC  REPORT  OCTOBER  20,1991 


INPOT  PAIRS  OF  DI.AI 


.050 

.035 

.100 

.063 

.150 

.080 

.200 

.100 

.300 

.120 

.400 

.140 

.500 

.160 

.700 

.170 

.800 

.200 

1.000 

.210 

2.000 

.270 

5.000 

.360 

10.000 

.400 

20.000 

.490 

50.000 

.590 

100.000 

.640 

200.000 

.700 

500.000 

.8001000.000 

.860 

INPUT 

.758 

PAIRS  OF  D(iiun)  ,P 
1.000  .642 

.992 

.543 

.962 

.460 

.880 

.390 

.721  .330 

.502 

.279 

.282 

.237 

.124 

.200 

.042  .170 

.012 

.144 

.000 

I 

Y(I) 

HO(I) 

HP(I) 

HEQ(l) 

AC(I) 

VUSEDC(I) 

PCC(I) 

1 

.000 

-1.500 

-1.500 

.000 

.000 

.000 

.000 

2 

1.000 

-1.400 

-1.500 

.000 

.000 

.000 

.000 

3 

2.000 

-1.300 

-1.500 

.000 

.000 

.000 

.000 

4 

3.000 

-1.200 

-1.500 

.000 

.000 

.000 

.000 

5 

4.000 

-1.100 

-1.500 

.000 

.000 

.000 

.000 

6 

5.000 

-1.000 

-1.500 

.000 

.000 

.000 

.000 

7 

6.000 

-.900 

-1.500 

.000 

.000 

.000 

.000 

8 

7.000 

-.800 

-1.500 

.000 

.000 

.000 

.000 

9 

8.000 

-.700 

-1.500 

.000 

.000 

.000 

.000 

10 

9.000 

-.600 

-1.500 

.000 

.000 

.000 

.000 

11 

10.000 

-.500 

-1.500 

.000 

.000 

.000 

.000 

12 

11.000 

-.400 

-1.500 

.000 

.000 

.000 

.000 

13 

12.000 

-.300 

-1.500 

.000 

.000 

.000 

.000 

14 

13.000 

-.200 

-1.500 

.000 

.000 

.000 

.000 

15 

14.000 

-.100 

-1.500 

.000 

.000 

.000 

.000 

16 

15.000 

.000 

-1.500 

.000 

.000 

.000 

.000 

17 

16.000 

.057 

-1.500 

.000 

.000 

.000 

.000 

18 

17.000 

.  100 

-1.500 

.000 

.000 

.000 

.000 

19 

18.000 

.  138 

-1.500 

.000 

.000 

.  000 

.000 

20 

19.000 

.  173 

-1.500 

.000 

.000 

.000 

.000 

21 

20.000 

.205 

-1.500 

.000 

.000 

.000 

.000 

22 

21.000 

.236 

-1.500 

-1.500 

.000 

.  000 

.000 

23 

22.000 

.266 

-1.500 

-1.414 

.000 

.000 

.000 

24 

23.000 

.295 

-1.500 

-1.364 

.000 

.000 

.000 

25 

24.000 

.  322 

-1.500 

-1.314 

.000 

.000 

.000 

26 

25.000 

.349 

-1.500 

-1.264 

.000 

.000 

.000 

27 

26.000 

.375 

-1.500 

-1.214 

.000 

.000 

.000 

28 

27.000 

.400 

-1.500 

-1.164 

.000 

.000 

.000 

29 

28.000 

.425 

-1.500 

-1.114 

.000 

.000 

.000 

30 

29.000 

.449 

-1.500 

-1.064 

.000 

.000 

.000 

31 

30.000 

.473 

-1.500 

-1.014 

.000 

.000 

.000 

32 

31.000 

.496 

-1.500 

-.964 

.000 

.000 

.000 

33 

32.000 

.519 

-1.500 

-.914 

.000 

.000 

.000 

34 

33.000 

.541 

-1.500 

-.864 

.000 

.000 

.000 

35 

34.000 

.563 

-1.500 

-.814 

.000 

.000 

.000 

36 

35.000 

.585 

-1.500 

-.764 

.000 

.000 

.000 

37 

36.000 

.606 

-1.500 

-.714 

.000 

.000 

.000 

38 

37.000 

.628 

-1.500 

—  .664 

.000 

.000 

.000 

39 

38.000 

.649 

-1.500 

-.614 

.000 

.000 

.000 
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40 

39.000 

.  669 

-1.500 

-.564 

.000 

.000 

.000 

41 

40.000 

.689 

-1.500 

-.514 

.000 

.000 

.000 

42 

41.000 

.710 

-1.500 

-.464 

.000 

.000 

.000 

43 

42.000 

.730 

-1.500 

-.414 

.000 

.000 

.000 

44 

43.000 

.749 

-1.500 

-.364 

.000 

.000 

.000 

45 

44.000 

.769 

-1.500 

-.314 

.000 

.000 

.000 

46 

45.000 

.788 

-1.500 

-.264 

.000 

.000 

.000 

47 

46.000 

.807 

-1.500 

-.214 

.000 

.000 

.000 

48 

47.000 

.826 

-1.500 

-.164 

.000 

.000 

.000 

49 

48.000 

.845 

-1.500 

-.114 

.000 

.000 

.000 

50 

49.000 

.863 

-1.500 

-.064 

.000 

.000 

.000 

51 

50.000 

.882 

-1.500 

-.014 

.156 

6.346 

.090 

52 

51.000 

.900 

-1.500 

.032 

.155 

6.546 

.093 

53 

52.000 

.918 

-1.500 

.071 

,155 

6.746 

.096 

54 

53.000 

.936 

-1.500 

.107 

.155 

6.946 

.099 

55 

54.000 

.954 

-1.500 

.142 

.154 

7.146 

.102 

56 

55.000 

.972 

-1.500 

.176 

.154 

7.346 

.104 

57 

56.000 

.989 

-1-500 

.208 

.154 

7.546 

.107 

58 

57.000 

1.007 

-1.500 

.239 

.153 

7.746 

.110 

59 

58.000 

1.024 

-1-500 

.270 

.153 

7.946 

.113 

60 

59.000 

1.041 

-1.429 

.300 

.153 

8.146 

.116 

61 

60.000 

1.058 

-1-329 

.329 

.152 

8.346 

.119 

62 

61.000 

1.075 

-1.229 

.358 

.152 

8.546 

.121 

63 

62.000 

1.092 

-1.129 

.386 

.152 

8.746 

.124 

64 

63.000 

1.109 

-1.029 

.414 

.151 

8.946 

.127 

65 

64.000 

1.125 

-.929 

.441 

.151 

9.146 

.130 

66 

65,000 

1.142 

-.829 

.468 

.151 

9.346 

.133 

67 

66.000 

1.158 

-.729 

.495 

.151 

9.546 

.136 

68 

67.000 

1-175 

-.629 

.521 
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.000 

.000 

739 

738.000 

7.489 

.000 

.000 

.000 

.000 

.000 

740 

739.000 

7.496 

.000 

.000 

.000 

.  000 

.000 

741 

740.000 

7.503 

.000 

.000 

.000 

.000 

.000 

742 

741.000 

7.511 

.000 

.000 

.000 

.000 

.000 

743 

742.000 

7.518 

.000 

.000 

.000 

.000 

.000 

744 

743.000 

7.525 

.000 

.000 

.000 

.000 

.000 

745 

744.000 

7.532 

.000 

.000 

.000 

.000 

.000 

746 

745.000 

7.539 

.000 

.000 

.000 

.000 

.000 

747 

746.000 

7.546 

.000 

.000 

.000 

.000 

.000 

748 

747.000 

7.553 

.000 

.000 

.000 

.000 

.000 

749 

748.000 

7.561 

.000 

.000 

.000 

.000 

.000 

750 

749.000 

7.568 

.000 

.000 

.000 

.000 

.000 

751 

750.000 

7.575 

.000 

.000 

.000 

.000 

.000 

752 

751.000 

7.582 

.000 

.000 

.000 

.000 

.000 

753 

752.000 

7.589 

.000 

.000 

.000 

.000 

.000 

754 

753.000 

7.596 

.000 

.000 

.000 

.000 

.000 

755 

754.000 

7.603 

.000 

.000 

.000 

.000 

.000 

756 

755.000 

7.610 

.000 

.000 

.000 

.000 

.000 

757 

756.000 

7.617 

.000 

.000 

.000 

.  000 

.000 

758 

757.000 

7.624 

.000 

.000 

.000 

.000 

.000 

759 

758.000 

7.632 

.000 

.000 

.000 

.000 

.000 
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760 

759.000 

7.639 

.000 

761 

760.000 

7.646 

.000 

762 

761.000 

7.653 

.000 

763 

762.000 

7.660 

.000 

764 

763.000 

7 . 667 

.000 

765 

764.000 

7.674 

.000 

766 

765.000 

7-681 

.000 

767 

766.000 

7.688 

.000 

768 

767.000 

7.695 

.000 

769 

768.000 

7-702 

.000 

770 

769.000 

7.709 

.000 

771 

770.000 

7.716 

.000 

772 

771.000 

7.723 

.000 

773 

772 . 000 

7.731 

.000 

774 

773.000 

7.738 

.000 

775 

774.000 

7.745 

.000 

776 

775.000 

7.752 

.000 

777 

776.000 

7.759 

.000 

778 

777.000 

7.766 

.000 

779 

778.000 

7.773 

.000 

780 

779.000 

7.780 

.000 

I 

DYEQ(I) 

VGENV(l) 

VUSED(I 

1 

20.000 

71.105 

55.379 

2 

21.000 

68.445 

276.944 

3 

20.500 

69.787 

150.520 

4 

20.000 

71.105 

55.387 

5 

20.250 

70.445 

67.165 

6 

20.500 

69.787 

150.504 

7 

20.375 

70.116 

106.640 

8 

20.250 

70.445 

67.168 

9 

20.313 

70.280 

83.041 

10 

20.281 

70.363 

72.685 

11 

20.250 

70.445 

67.166 

12 

20.266 

70.404 

69.207 

13 

20.281 

70.363 

72.684 

14 

20.273 

70.383 

70.681 

15 

20.266 

70.404 

69.208 

16 

20.270 

70.394 

69.862 

17 

20.273 

70.383 

70.681 

18 

20.271 

70.388 

70.251 

19 

20.272 

70.386 

70.458 

20 

20.272 

70.387 

70.353 

21 

20.272 

70.386 

70.408 

22 

20.272 

70  387 

70.382 

23 

20.272 

70.387 

70.392 

24 

20.272 

70.387 

70.387 

25 

20.272 

70.387 

70.382 

26 

20.272 

70.387 

70.386 

27 

20.272 

70-387 

70.387 

28 

20.272 

70.387 

70.386 

29 

20.272 

70.387 

70.386 

30 

20.272 

70.387 

70-387 

31 

20.272 

70.387 

70.386 

32 

20.272 

70.387 

70.387 

33 

20.272 

70.387 

70.387 

34 

20.272 

70.387 

70.387 

35 

20.272 

70.387 

70.387 

Figure  A3.  (Sheet  14  of  15) 

A26 


.000  .000  .000  .000 

.000  .000  .000  .000 

.000  .000  .000  .000 

.000  .000  .000  .000 

.000  .000  .000  .000 

.000  .000  .000  .000 

.000  .000  .000  .000 

.000  .000  .000  .000 

.000  .000  .000  .000 

.000  .000  .000  .000 

.000  .000  .000  .000 

.000  .000  .000  .000 

.000  .000  .000  .000 

.000  .000  .000  .000 

.000  .000  .000  .000 

.000  .000  .000  .000 

.000  .000  .000  .000 

.000  .000  .000  .000 

.000  .000  .000  .000 

.000  .000  .000  .000 

000  .000  .000  .000 


PSV ( I ) 
.508 
.489 
.498 
.508 
.503 
.498 
.501 
.503 
.502 
.503 
.503 
.503 
.503 
.503 
.503 
.503 
.503 
.503 
.503 
.503 
.503 
.503 
.503 
.503 
.503 
.503 
.503 
.503 
.503 
.503 
.503 
.503 
.503 
.503 
.503 
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36 

20.272 

70.387 

70.387 

.503 

37 

20.272 

70.387 

70.387 

.503 

38 

20.272 

70.387 

70.387 

.503 

39 

20.272 

70.387 

70.387 

.503 

40 

20.272 

70.387 

70.387 

.503 

SOLUTION  REACHED,  NON-INTERSECTING  PROFILES 
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EXAMPLE  4,  CERC  REPORT  OCTOBER  23,1991  J 

^T-ryt'S  ^s«  ,sP  ^3TA«.  ^h«.h 

1  780  1.0  0.10  0.10  0.10  600.0  0.10  -1.5  6.0  0.2 


9  • 

NUMBER  OF  INPUT  (Y,HO) 

POINTS 

ON  ORIGINAL  PROFILE 

0.0 

-1.5 

15.0 

0.0 

90.0 

1.0 

140.0 

,  n  1  »•< 

t'n  t  W*'*’ 

200.0 

3.0 

300.0 

4.0 

415.0 

5.0 

540.0 

770.0 

^JSas 

7.0 

^STc,  Chl»1t. 

;  riKL/,sr« 

in  Thif 

1 

3  2. 

64  1 

0.10 

JSED, 

IMAX,  & 

XMU,SIG 

in  phi  units 

r - - 

.  — - 

-  -  - 

—  —  —  —  - 

—  —  - 

—  -  — 

‘  —  —  —  - 

'  oi,Ar 

•  0.05 

0.035 

0.10 

0.063 

0.15 

0.08 

0.20 

0 . 100  1  ( |i|  ^irt^ 

1  0.30 

0.120 

0.40 

0.140 

0.50 

0.160 

0.70 

0.170  'Tfc.iinput 

0*36  '  Same  ^ 

1  0.80 

0.200 

1.00 

0.210 

2.0 

0.27 

5.0 

1  10.0 

0.40 

20.0 

0.49 

50.0 

0.59 

100.0 

;p«  R-«» 

1  200.0 

0.70 

500.0 

0.80 

1000.0 

0.86 

C-  —  —  _ 

.  - - - 

—  — - 

—  —  -  -  — 

• 

— “  —  —  — 

—  -  —  — 

0.6 

xnifttK, 

1.0 

pye^i 

0.35 

0.5 

H 

O 

o 

• 

o 

3  pK-irs 

40  100.0 


A/oic.  ■  Lln&i  flfeovc  /-/ovc.  "feecrv 

"ho  Fft-cW'+cht  AnyV3'^«aiW'n. 
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EXAMPLE  4,  CERC  REPORT  OCTOBER  23,1991 


INPOT  PAIRS  OF  DI,AI 


.050 

.035 

.100 

.063 

.150 

.080 

.200 

.100 

.300 

.120 

.400 

.140 

.500 

.160 

.700 

.170 

.800 

.200 

1.000 

.210 

2.000 

.270 

5.000 

.360 

10.000 

200.000 

.400 

.700 

20.000 

500.000 

.490  50.000 

.8001000.000 

.590 

860 

100.000 

.640 

INPUT  PAIRS  OF  D(inm),P 

.600  1.000  .350 

.500 

.100 

.000 

I 

Y(I) 

HO(l) 

HP(I) 

HEQ(I) 

AC(I) 

VUSEDC(I) 

PCC(I) 

1 

.000 

-1.500 

-1.500 

.000 

.000 

.000 

.000 

2 

1.000 

-1.400 

-1.500 

.000 

.000 

.000 

.000 

3 

2.000 

-1.300 

-1.500 

.000 

.000 

.000 

.000 

4 

3.000 

-1.200 

-1.500 

.000 

.000 

.000 

.000 

5 

4.000 

-1.100 

-1.500 

.000 

.000 

.000 

.000 

6 

5.000 

-1.000 

-1-500 

.000 

.000 

.000 

.000 

7 

6.000 

-.900 

-1.500 

.000 

.000 

.000 

.000 

8 

7.000 

-.800 

-1-500 

.000 

.000 

.000 

.000 

9 

8.030 

-.700 

-1.500 

.000 

.000 

.000 

.000 

10 

9.000 

-.600 

-1.500 

.000 

.000 

.000 

.000 

11 

10.000 

-.500 

-1.500 

.000 

.000 

>00 

.000 

12 

11.000 

-.400 

-1.500 

.000 

.000 

.000 

.000 

13 

12.000 

-.300 

-1.500 

.000 

.000 

.000 

.000 

14 

13.000 

-.200 

-1.500 

.000 

.000 

.000 

.000 

15 

14.000 

-.100 

-1.500 

.000 

.000 

.000 

.000 

16 

15.000 

.000 

-1.500 

.000 

.000 

.000 

.000 

17 

16.000 

.013 

-1.500 

.000 

.000 

.000 

.000 

18 

17.000 

.027 

-1.500 

.000 

.000 

.000 

.000 

19 

18.000 

.040 

-1.500 

.000 

.000 

.000 

.000 

20 

19.000 

.053 

-1.500 

.000 

.000 

.000 

.000 

21 

20.000 

.067 

-1.500 

.000 

.000 

.000 

.000 

22 

21.000 

.080 

-1.500 

.000 

.000 

.000 

.000 

23 

22.000 

.093 

-1.500 

.000 

.000 

.000 

.000 

24 

23.000 

.107 

-1.500 

.000 

.000 

.000 

.000 

25 

24.000 

.120 

-1.500 

.000 

.000 

.000 

.000 

26 

25.000 

.133 

-1.500 

.000 

.000 

.000 

.000 

27 

26.000 

.147 

-1.500 

.000 

.000 

.000 

.000 

28 

27.000 

.160 

-1.500 

.000 

.000 

.000 

.000 

29 

28.000 

.173 

-1.500 

.000 

.000 

.000 

.000 

30 

29.000 

.187 

-1.500 

.000 

.000 

.000 

.000 

31 

30.000 

.200 

-1.500 

.000 

.000 

.000 

.000 

32 

31.000 

.213 

-1.500 

.000 

.000 

.000 

.000 

33 

32.000 

.227 

-1.500 

.000 

.000 

.000 

.000 

34 

33.000 

.240 

-1.500 

.000 

.000 

.000 

.000 

35 

34.000 

.253 

-1.500 

.000 

.000 

.000 

.000 

36 

35.000 

.267 

-1.500 

.000 

.000 

.000 

.000 

37 

36.000 

.280 

-1.500 

.000 

.000 

.000 

.000 

38 

37.000 

.293 

-1.500 

.000 

.000 

.000 

.000 

39 

38.000 

.307 

-1.500 

.000 

.000 

.000 

.000 

40 

39.000 

.320 

-1.500 

.000 

.000 

.000 

.000 

41 

40.000 

.333 

-1.500 

.000 

.000 

.000 

.000 
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42 

41.000 

.347 

-1.500 

.000 

.000 

.000 

.000 

43 

42.000 

.360 

-1.500 

.000 

.000 

.000 

.000 

44 

43.000 

.373 

-1.500 

.000 

.000 

.000 

.000 

45 

44.000 

.387 

-1.500 

.000 

.000 

.000 

.000 

46 

45.000 

.400 

-1.500 

.000 

.000 

.000 

.000 

47 

46.000 

.413 

-1.500 

.000 

.000 

.000 

.000 

48 

47.000 

.427 

-1.500 

.000 

.000 

.000 

.000 

49 

48.000 

.440 

-1.500 

.000 

.000 

.000 

.000 

50 

49.000 

.453 

-1.500 

.000 

.000 

.000 

.000 

51 

50.000 

.467 

-1.500 

.000 

.000 

.000 

.000 

52 

51.000 

.480 

-1.500 

.000 

.000 

.000 

.000 

53 

52.000 

.493 

-1.500 

.000 

.000 

.000 

.000 

54 

53.000 

.507 

-1.500 

.000 

.000 

.000 

.000 

55 

54.000 

.520 

-1.500 

.000 

.000 

.000 

.000 

56 

55.000 

.533 

-1.500 

.000 

.000 

.000 

.000 

57 

56.000 

.547 

-1.500 

.000 

.000 

.000 

.000 

58 

57.000 

.560 

-1.500 

.000 

.000 

.000 

.000 

59 

58.000 

.573 

-1.500 

.000 

.000 

.000 

.000 

60 

59.000 

.587 

-1.500 

.000 

.000 

.000 

.000 

61 

60.000 

.600 

-1.500 

.000 

.000 

.000 

.000 

62 

61.000 

.613 

-1.500 

.000 

.000 

.000 

.000 

63 

62.000 

.627 

-1.500 

.000 

.000 

.000 

.000 

64 

63.000 

.640 

-1.500 

.000 

.000 

.000 

.000 

65 

64.000 

.653 

-1.500 

.000 

.000 

.000 

.000 

66 

65.000 

.667 

-1.500 

.000 

.000 

.000 

.000 

67 

66.000 

.680 

-1.500 

.000 

.000 

.000 

.000 

68 

67.000 

.693 

-1.500 

.000 

.000 

.000 

.000 

69 

68.000 

.707 

-1.500 

.000 

.000 

.000 

.000 

70 

69.000 

.720 

-1.500 

.000 

.000 

.000 

.000 

71 

70.000 

.733 

-1.500 

.000 

.000 

.000 

.000 

72 

71.000 

.747 

-1.500 

.000 

.000 

.000 

.000 

73 

72.000 

.760 

•1.500 

.000 

.000 

.000 

.000 

74 

73.000 

.773 

-1.500 

.000 

.000 

.000 

.000 

75 

74.000 

.787 

-1.500 

.000 

.000 

.000 

.000 

76 

75.000 

.800 

-1.500 

.000 

.000 

.000 

.000 

77 

76.000 

.813 

-1.500 

.000 

.000 

.000 

.000 

78 

77.000 

.827 

-1.500 

.000 

.000 

.000 

.000 

79 

78.000 

.840 

-1.500 

.000 

.000 

.000 

.000 

80 

79.000 

.853 

-1.500 

.000 

.000 

.000 

.000 

81 

80.000 

.867 

-1.500 

.000 

.000 

.000 

.000 

82 

81.000 

.880 

-1.500 

.000 

.000 

.000 

.000 

83 

82.000 

.893 

-1.500 

.000 

.000 

.000 

.000 

84 

83.000 

.907 

-1.500 

.000 

.000 

.000 

.000 

85 

84.000 

.920 

-1.500 

.000 

.000 

.000 

.000 

86 

85.000 

.933 

-1.500 

.000 

.000 

.000 

.000 

87 

86.000 

.947 

-1.500 

.000 

.000 

.000 

.000 

88 

87.000 

.960 

-1.500 

.000 

.000 

.000 

.000 

89 

88.000 

.973 

-1.500 

.000 

.000 

.000 

.000 

90 

89.000 

.987 

-1.500 

.000 

.000 

.000 

.000 

91 

90.000 

1.000 

-1.500 

.000 

.000 

.000 

.000 

92 

91.000 

1.020 

-1.500 

.000 

.000 

.000 

.000 

93 

92.000 

1.040 

-1.500 

.000 

.000 

.000 

.000 

94 

93.000 

1.060 

-1.500 

.000 

.000 

.000 

.000 

95 

94.000 

1.080 

-1.500 

.000 

.000 

.000 

.000 

96 

95.000 

1.100 

-1.500 

.000 

.000 

.000 

.000 

97 

96.000 

1.120 

-1.500 

.000 

.000 

.000 

.000 

98 

97.000 

1.140 

-1.500 

.000 

.000 

.000 

.000 

99 

98.000 

1.160 

-1.500 

.000 

.000 

.000 

.000 

100 

99.000 

1.180 

-1.500 

.000 

.000 

.000 

.000 

101 

100.000 

1.200 

-1.500 

.000 

.000 

.000 

.000 
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102 

101.000 

1.220 

-1.500 

103 

102.000 

1.240 

-1.500 

104 

103.000 

1.260 

-1.500 

105 

104.000 

1.280 

-1.500 

106 

105.000 

1.300 

-1.500 

107 

106.000 

1.320 

-1.500 

108 

107.000 

1.340 

-1.500 

109 

108.000 

1.360 

-1.500 

110 

109.000 

1.380 

-1.500 

111 

110.000 

1.400 

-1.500 

112 

111.000 

1.420 

-1.500 

113 

112.000 

1.440 

-1.500 

114 

113.000 

1.460 

-1.500 

115 

114.000 

1.480 

-1.500 

116 

115,000 

1.500 

-1.500 

117 

116.000 

1.520 

-1.500 

118 

117.000 

1.540 

-1.500 

119 

118.000 

1.560 

-1.500 

120 

119.000 

1.580 

-1.500 

121 

120.000 

1.600 

-1.500 

122 

121.000 

1.620 

-1.500 

123 

122.000 

1.640 

-1.500 

124 

123.000 

1.660 

-1.500 

125 

124.000 

1.680 

-1.500 

126 

125.000 

1.700 

-1.500 

127 

126.000 

1.720 

-1.500 

128 

127.000 

1.740 

-1.500 

129 

128.000 

1.760 

-1.500 

130 

129.000 

1.780 

-1.500 

131 

130.000 

1.800 

-1.500 

132 

131.000 

1.820 

-1.500 

133 

132.000 

1.840 

-1,500 

134 

133.000 

1.860 

-1.500 

135 

134.000 

1.880 

-1.500 

136 

135.000 

1.900 

-1,500 

137 

136.000 

1.920 

-1.500 

138 

137.000 

1.940 

-1,500 

139 

138.000 

1.960 

-1.500 

140 

139.000 

1.980 

-1.500 

141 

140.000 

2.000 

-1.500 

142 

141.000 

2.017 

-1,500 

143 

142.000 

2.033 

-1.500 

144 

143.000 

2.050 

-1.500 

145 

144.000 

2.067 

-1,500 

146 

145.000 

2.083 

-1.500 

147 

146.000 

2.100 

-1,500 

148 

147,000 

2.117 

-1.500 

149 

148.000 

2,133 

-1.500 

150 

149.000 

2.150 

-1.500 

151 

150.000 

2.167 

-1.500 

152 

151.000 

2.183 

-1.500 

153 

152.000 

2.200 

-1.500 

154 

153.000 

2.217 

-1.500 

155 

154.000 

2.233 

-1.500 

156 

155.000 

2.250 

-1.500 

157 

156.000 

2.267 

-1.500 

158 

157.000 

2.283 

-1.500 

159 

158.000 

2.300 

-1.500 

160 

159.000 

2.317 

-1.500 

161 

160.000 

2.333 

-1.500 

Figure  A5.  (Sheet  3  of  15) 


.000  .000  .000  .000 

.000  .000  .000  .000 

.000  .000  .000  .000 

.000  .000  .000  .000 

.000  .000  .000  .000 

.000  .000  .000  .000 

.000  .000  .000  .000 

.  000  . 000  . 000  . 000 

.000  .000  .000  .000 

.000  .000  .000  .000 

.000  .000  .000  .000 

.000  .000  .000  .000 

.000  .000  .000  .000 

.000  .000  .000  .000 

.000  .000  .000  .000 

.000  .000  .000  .000 

.000  .000  .000  .000 

.000  .000  .000  .000 

.000  .000  .000  .000 

.000  .000  .000  .000 

.000  .000  .000  .000 

.000  .000  .000  .000 

.000  .000  .000  .000 

.000  .000  .000  .000 

.000  .000  .000  .000 

.000  .000  .000  .000 

.000  .000  .000  .000 

.000  .000  .000  .000 

.000  .000  .000  .000 

.000  .000  .000  .000 

.000  .000  .000  .000 

.000  .000  .000  .000 

.000  .000  .000  .000 

.000  .000  .000  .000 

.000  .000  .000  .000 

.000  .000  .000  .000 

.000  .000  .000  .000 

.000  .000  .000  .000 

.000  .000  .000  .000 

-1.500  .000  .000  .000 

-1.370  ,000  .000  .000 

-1.270  .000  .000  .000 

-1.170  .000  ,000  .000 

-1.070  .000  .000  .000 

-.970  .000  .000  .000 

-.870  .000  .000  .000 

-.770  .000  .000  .000 

-.670  .000  .000  .000 

-.570  .000  .000  .000 

-.470  .000  .000  .000 

-.370  .000  .000  .000 

-.270  .000  .000  .000 

-.170  .000  .000  .000 

-.070  .164  3.060  .021 

.024  .163  3.260  .022 

.089  .163  3.460  .023 

.145  .163  3.660  .025 

.197  .163  3.860  .026 

.245  .163  4.060  .027 

.291  .163  4.260  .029 


Appondix  A  Listing  of  Program  EQPR.FOR 


A31 


162 

161.000 

2.350 

-1.500 

.335 

.163 

4.460 

163 

162.000 

2.367 

-1.500 

.377 

.163 

4.660 

164 

163.000 

2.383 

-1.500 

.418 

.163 

4.860 

165 

164.000 

2.400 

-1.500 

.458 

.163 

5.060 

166 

165.000 

2.417 

-1.500 

.497 

.163 

5.260 

167 

166.000 

2.433 

-1.500 

.535 

.162 

5.460 

168 

167.000 

2.450 

-1.500 

.572 

.162 

5.660 

169 

168.000 

2.467 

-1.500 

.608 

.162 

5.860 

170 

169.000 

2.483 

-1.500 

.643 

.162 

6.060 

171 

170.000 

2.500 

-1.500 

.678 

.162 

6.260 

172 

171.000 

2.517 

-1.500 

.712 

.162 

6.460 

173 

172.000 

2.533 

-1.500 

.746 

.162 

6.660 

174 

173.000 

2.550 

-1.500 

.779 

.162 

6.860 

175 

174.000 

2.567 

-1.500 

.812 

.162 

7.060 

176 

175.000 

2.583 

-1.500 

.844 

.162 

7.260 

177 

176.000 

2.600 

-1.500 

.876 

.161 

7.460 

178 

177.000 

2.617 

-1.500 

.908 

.161 

7.660 

179 

178.000 

2.633 

-1.500 

.939 

.161 

7.860 

180 

179.000 

2.650 

-1.500 

.969 

.161 

8.060 

181 

180.000 

2.667 

-1.500 

1.000 

.161 

8.260 

182 

181.000 

2.683 

-1.500 

1.030 

.161 

8.460 

183 

182.000 

2.700 

-1.500 

1.059 

.161 

8.660 

184 

183.000 

2.717 

-1.500 

1.089 

.161 

8.860 

185 

184.000 

2.733 

-1.500 

1.118 

.161 

9.060 

186 

185.000 

2.750 

-1.500 

1.146 

.161 

9.260 

187 

186.000 

2.767 

-1.500 

1.175 

.161 

9.460 

188 

187.000 

2.783 

-1.500 

1.203 

.160 

9.660 

189 

188.000 

2.800 

-1.406 

1.231 

.160 

9.860 

190 

189.000 

2.817 

-1.306 

1.259 

.160 

10.060 

191 

190.000 

2.833 

-1.206 

1.286 

.160 

10.260 

192 

191.000 

2.850 

-1.106 

1.313 

.160 

10.460 

193 

192.000 

2.867 

-1.006 

1.343. 

.160 

10.660 

194 

193.000 

2.883 

-.906 

1.367 

.160 

10.860 

195 

194.000 

2.900 

-.806 

1.394 

.160 

11.060 

196 

195.000 

2.917 

-.706 

1.420 

.160 

11.260 

197 

196.000 

2.933 

-.606 

1.447 

.160 

11.460 

198 

197.000 

2.950 

-.506 

1.473 

.159 

11.660 

199 

198.000 

2.967 

-.406 

1.498 

.159 

11.860 

200 

199.000 

2.983 

-.306 

1.524 

.159 

12.060 

201 

200.000 

3.000 

-.206 

1.550 

.159 

12.260 

202 

201.000 

3.010 

-.106 

1.575 

.159 

12.460 

203 

202.000 

3.020 

-.006 

1.600 

.159 

12.660 

204 

203.000 

3.030 

.094 

1.625 

.159 

12.860 

205 

204.000 

3.040 

.194 

1.650 

.159 

13.060 

206 

205.000 

3.050 

.294 

1.674 

.159 

13.260 

207 

206.000 

3.060 

.394 

1-699 

.159 

13.460 

208 

207.000 

3.070 

.494 

1.723 

.159 

13.660 

209 

208.000 

3.080 

.594 

1.748 

.158 

13.860 

210 

209.000 

3.090 

.694 

1.772 

.158 

14.060 

211 

210.000 

3.100 

.794 

1.796 

.158 

14.260 

212 

211.000 

3.110 

.894 

1.819 

.158 

14.460 

213 

212.000 

3.120 

.994 

1.843 

.158 

14.660 

214 

213.000 

3.130 

1.094 

1.867 

.158 

14.860 

215 

214.000 

3.140 

1.194 

1.890 

.158 

15.060 

216 

215.000 

3.150 

1.294 

1.913 

.158 

15.260 

217 

216.000 

3.160 

1.394 

1.936 

.158 

15.460 

218 

217.000 

3.170 

1.494 

1.959 

.158 

15.660 

219 

218.000 

3.180 

1.594 

1.982 

.158 

15.860 

220 

219.000 

3.190 

1.694 

2.005 

.157 

16.060 

221 

220.000 

3.200 

1.794 

2.028 

.157 

16.260 

Figure  A5.  (Sheet  4  of  15) 


A32 


.030 

.031 

.033 

.034 

.036 

.037 

.038 

.040 

.041 

.042 

.044 

.045 

.046 

.048 

.049 

.050 

.052 

.053 

.054 

.056 

.057 

.058 

.060 

.061 

.063 

.064 

.065 

.067 

.068 

.069 

.071 

.072 

.073 

.075 

.076 

.077 

.079 

.080 

.081 

.083 

.084 

.086 

.087 

.088 

.090 

.091 

.092 

.094 

.095 

.096 

.098 

.099 

.100 

.102 

.103 

.104 

.106 

.107 

.108 

.110 


Appendix  A  Listing  of  Program  EQPR.FOR 


222 

221.000 

3.210 

1.894 

223 

222.000 

3.220 

1.994 

224 

223.000 

3.230 

2.094 

225 

224.000 

3.240 

2.194 

226 

225.000 

3.250 

2.294 

227 

226.000 

3.260 

2.394 

228 

227.000 

3.270 

2.494 

229 

228.000 

3.280 

2.594 

230 

229.000 

3.290 

2.694 

231 

230.000 

3.300 

2.794 

232 

231.000 

3.310 

2.894 

233 

232.000 

3.320 

2.994 

234 

233.000 

3.330 

3.094 

235 

234.000 

3.340 

3.194 

236 

235.000 

3.350 

3.294 

237 

236.000 

3.360 

3.394 

238 

237.000 

3.370 

.000 

239 

238.000 

3.380 

.000 

240 

239.000 

3.390 

.000 

241 

240.000 

3.400 

.000 

242 

241.000 

3.410 

.000 

243 

242.000 

3.420 

.000 

244 

243.000 

3.430 

.000 

245 

244.000 

3.440 

.000 

246 

245.000 

3.450 

.000 

247 

246.000 

3.460 

.000 

248 

247.000 

3.470 

.000 

249 

248.000 

3.480 

.000 

250 

249.000 

3.490 

.000 

251 

250.000 

3.500 

.000 

252 

251.000 

3.510 

.000 

253 

252.000 

3.520 

.000 

254 

253.000 

3.530 

.000 

255 

254.000 

3.540 

.000 

256 

255.000 

3.550 

.000 

257 

256.000 

3.560 

.000 

258 

257.000 

3.570 

.000 

259 

258.000 

3.580 

.000 

260 

259.000 

3.590 

.000 

261 

260.000 

3.600 

.000 

262 

261.000 

3.610 

.000 

263 

262.000 

3.620 

.000 

264 

263.000 

3.630 

.000 

265 

264.000 

3.640 

.000 

266 

265.000 

3.650 

.000 

267 

266.000 

3.660 

.000 

268 

267.000 

3.670 

.000 

269 

268.000 

3.680 

.000 

270 

269.000 

3.690 

.000 

271 

270.000 

3.700 

.000 

272 

271.000 

3.710 

.000 

273 

272.000 

3.720 

.000 

274 

273.000 

3.730 

.000 

275 

274.000 

3.740 

.000 

276 

275.000 

3.750 

.000 

277 

276.000 

3.760 

.000 

278 

277.000 

3.770 

.000 

279 

278.000 

3.780 

.000 

280 

279.000 

3.790 

.000 

281 

280.000 

3.800 

.000 

Figure  A5.  (Sheet  5  of  15) 

Appendix  A  Listing  of  Program  EQPR.FOR 


2.050 

.157 

16.460 

.111 

2.073 

.157 

16.660 

.113 

2.095 

.157 

16.860 

.114 

2.117 

.157 

17.101 

.116 

2.139 

.157 

17.216 

.116 

2.161 

.157 

17.409 

.118 

2.183 

.157 

17.680 

.119 

2.205 

.156 

18.030 

.122 

2.227 

.156 

18.457 

.125 

2.248 

.156 

18.963 

.128 

2.270 

.156 

19.547 

.132 

2.291 

.155 

20.210 

.137 

2.313 

.155 

20.952 

.142 

2.334 

.155 

21.772 

.147 

2.355 

.154 

22.672 

.153 

2.375 

.154 

23.435 

.158 

2.396 

.153 

24.414 

.165 

2.417 

.153 

25.383 

.171 

2.437 

.153 

26.341 

.178 

2.457 

.152 

27.289 

.184 

2.477 

.152 

28.227 

.191 

2.497 

.151 

29.155 

.197 

2.517 

.151 

30.072 

.203 

2.537 

.150 

30.980 

.209 

2.556 

.150 

31.879 

.215 

2.576 

.150 

32.768 

.221 

2.595 

.149 

33.647 

.227 

2.614 

.149 

34.517 

.233 

2.634 

.148 

35.378 

.239 

2.652 

.148 

36.230 

.245 

2.671 

.147 

37.074 

.250 

2.690 

.147 

37.908 

.256 

2.709 

.147 

38.733 

.262 

2.727 

.146 

39.550 

.267 

2.746 

.146 

40.359 

.273 

2.764 

.146 

41.159 

.278 

2.782 

.145 

41.951 

.283 

2.800 

.145 

42.735 

.289 

2.818 

.144 

43.511 

.294 

2.836 

.144 

44.279 

.299 

2.854 

.144 

45.039 

.304 

2.871 

.143 

45.791 

.309 

2.889 

.143 

46.536 

.314 

2.907 

.143 

47.273 

.319 

2.924 

.142 

48.003 

.324 

2.941 

.142 

48.726 

.329 

2.958 

.142 

49.441 

.334 

2.975 

.141 

50.149 

.339 

2.992 

.141 

50.850 

.343 

3.009 

.141 

51.544 

.348 

3.026 

.140 

52.231 

.353 

3.043 

.140 

52.912 

.357 

3.060 

.140 

53.585 

.362 

3.076 

.139 

54.252 

.366 

3.093 

.139 

54.913 

.371 

3.109 

.139 

55.567 

.375 

3.125 

.138 

56.215 

.380 

3.142 

.138 

56.856 

.384 

3.158 

.138 

57.491 

.388 

3.174 

.138 

58.120 

.393 

A33 


282 

281.000 

3.810 

.000 

283 

282.000 

3.820 

.000 

284 

283.000 

3.830 

.000 

285 

284 . 000 

3.840 

.000 

286 

285.000 

3.850 

.000 

287 

286.000 

3.860 

.000 

288 

287.000 

3.870 

.000 

289 

288.000 

3.880 

.000 

290 

289.000 

3.890 

.000 

291 

290.000 

3.900 

.000 

292 

291.000 

3.910 

.000 

293 

292.000 

3.920 

.000 

294 

293.000 

3.930 

.000 

295 

294.000 

3.940 

.000 

296 

295.000 

3.950 

.000 

297 

296.000 

3.960 

.000 

298 

297.000 

3.970 

.000 

299 

298.000 

3.980 

.000 

300 

299.000 

3.990 

.000 

301 

300.000 

4.000 

.000 

302 

301.000 

4.009 

.000 

303 

302.000 

4.017 

.000 

304 

303 . 000 

4.026 

.000 

305 

304.000 

4.035 

.000 

306 

305.000 

4.043 

.000 

307 

306.000 

4-052 

.000 

308 

307.000 

4.061 

.000 

309 

308.000 

4.070 

.000 

310 

309.000 

4.078 

.000 

311 

310.000 

4.087 

.000 

312 

311.000 

4.096 

.000 

313 

312.000 

4.104 

.000 

314 

313.000 

4.113 

.000 

315 

314.000 

4.122 

.000 

316 

315.000 

4.130 

.000 

317 

316.000 

4.139 

.000 

318 

317.000 

4.148 

.000 

319 

318.000 

4.157 

.000 

320 

319.000 

4.165 

.000 

321 

320.000 

4.174 

.000 

322 

321.000 

4.183 

.000 

323 

322.000 

4.191 

.000 

324 

323.000 

4.200 

.000 

325 

324.000 

4.209 

.000 

326 

325.000 

4.217 

,000 

327 

326.000 

4.226 

.000 

328 

327.000 

4.235 

.000 

329 

328.000 

4.243 

.000 

330 

329.000 

4.252 

.000 

331 

330.000 

4.261 

.000 

332 

331.000 

4.270 

.000 

333 

332.000 

4.278 

.000 

334 

333.000 

4.287 

,000 

335 

334.000 

4.296 

.000 

336 

335.000 

4.304 

.000 

337 

336.000 

4.313 

.000 

338 

337.000 

4.322 

.000 

339 

338.000 

4.330 

.000 

340 

339.000 

4.339 

.000 

341 

340.000 

4.348 

.000 

3-190 

.137 

58.743 

.397 

3.206 

.137 

59.361 

.401 

3.222 

.137 

59.972 

.405 

3.238 

.136 

60.577 

.409 

3.253 

.136 

61.177 

.413 

3.269 

.136 

61.771 

.417 

3.284 

.136 

62.359 

.421 

3.300 

.135 

62.942 

.425 

3.315 

.135 

63.519 

.429 

3.331 

.135 

64.091 

.433 

3.346 

.134 

64.658 

.437 

3.361 

.134 

65.219 

.440 

3.376 

.134 

65.776 

.444 

3.391 

.134 

66.327 

.448 

3.406 

.133 

66.873 

.452 

3.421 

.133 

67.414 

.455 

3.436 

.133 

67.950 

.459 

3.451 

.133 

68.482 

.463 

3.466 

.132 

69.008 

.466 

3.481 

.132 

69.530 

.470 

3.495 

.132 

70.047 

.473 

3.510 

.132 

70.557 

.477 

3.524 

.131 

71.062 

.480 

3.539 

.131 

71.561 

.483 

3.553 

.131 

72.054 

.487 

3.567 

.131 

72.542 

.490 

3.582 

.130 

73.024 

.493 

3.596 

,130 

73.500 

.496 

3.610 

.130 

73.971 

.500 

3.624 

.130 

74.436 

.503 

3.638 

.129 

74.896 

.506 

3.652 

.129 

75.351 

.509 

3.666 

.128 

75.801 

.512 

3.680 

.128 

76.245 

.515 

3.694 

.128 

76.684 

.518 

3.707 

.127 

77.119 

.521 

3.721 

.127 

77.548 

.524 

3.734 

.126 

77.973 

.527 

3.748 

.126 

78.393 

.529 

3.761 

.126 

78.808 

.532 

3.774 

.125 

79.218 

.535 

3.788 

.125 

79.624 

.538 

3.801 

.125 

80.026 

.541 

3.814 

.124 

80.423 

.543 

3.827 

.124 

80.816 

.546 

3.840 

.124 

81.204 

.548 

3.853 

.123 

81.589 

.551 

3.865 

.123 

81.969 

.554 

3.878 

.122 

82.345 

.556 

3.891 

.122 

82.717 

.559 

3,903 

.122 

83.086 

.561 

3.916 

.121 

83.450 

.564 

3.928 

.121 

83.810 

.566 

3.941 

.121 

84.167 

.568 

3.953 

.121 

84.520 

.571 

3.965 

.120 

84.870 

.573 

3.978 

.120 

85.216 

.576 

3.990 

.120 

85,558 

.578 

4.002 

.119 

85.897 

.580 

4.014 

.119 

86.233 

.582 

Figure  A5.  (Sheet  6  of  15) 


A34 


Appendix  A  Listing  of  Program  EQPR.FOR 


342 

341.000 

4.357 

.000 

4.026 

.119 

86.565 

.585 

343 

342.000 

4.365 

.000 

4.038 

.118 

86.894 

.587 

344 

343.000 

4.374 

.000 

4.050 

.118 

87.220 

.589 

345 

344.000 

4.383 

.000 

4.062 

.118 

87.542 

.591 

346 

345.000 

4.391 

.000 

4.073 

.117 

87.862 

.593 

347 

346.000 

4.400 

.000 

4.085 

.117 

88.178 

.596 

348 

347.000 

4.409 

.000 

4.097 

.117 

88.491 

.598 

349 

348.000 

4.417 

.000 

4.108 

.117 

88.802 

.600 

350 

349.000 

4.426 

.000 

4.120 

.116 

89.109 

.602 

351 

350.000 

4.435 

.000 

4.132 

.116 

89.414 

.604 

352 

351.000 

4.443 

.000 

4.143 

.116 

89.716 

.606 

353 

352.000 

4.452 

.000 

4.154 

.116 

90.015 

.608 

354 

353.000 

4.461 

.000 

4.166 

.115 

90.311 

.610 

355 

354.000 

4.470 

.000 

4.177 

.115 

90.605 

.612 

356 

355.000 

4.478 

.000 

4.188 

.115 

90.896 

.614 

357 

356.000 

4.487 

.000 

4.200 

.114 

91.185 

.616 

358 

357.000 

4.496 

.000 

4.211 

.114 

91.471 

.618 

359 

358.000 

4.504 

.000 

4.222 

.114 

91.754 

.620 

360 

359.000 

4.513 

.000 

4.233 

.114 

92.035 

.622 

361 

360.000 

4.522 

.000 

4.244 

.113 

92.314 

.623 

362 

361.000 

4.530 

.000 

4.255 

.113 

92.591 

.625 

363 

362.000 

4.539 

.000 

4.266 

.113 

92.865 

.627 

364 

363.000 

4.548 

.000 

4.277 

.113 

93.137 

.629 

365 

364.000 

4.557 

.000 

4.288 

.112 

93.407 

.631 

366 

365.000 

4.565 

.000 

4.299 

.112 

93.674 

.633 

367 

366.000 

4.574 

.000 

4.309 

.112 

93.940 

.634 

368 

367.000 

4.583 

.000 

4.320 

.112 

94.203 

.636 

369 

368.000 

4.591 

.000 

4.331 

.112 

94.465 

.638 

370 

369.000 

•«.600 

.000 

4.342 

.111 

94.724 

.640 

371 

370.000 

4.609 

.000 

4.352 

.111 

94.982 

.642 

372 

371.000 

4.617 

.000 

4.363 

.111 

95.237 

.643 

373 

372.000 

4.626 

.000 

4.373 

.111 

95.491 

.645 

374 

373.000 

4.635 

.000 

4.384 

.110 

95.743 

.647 

375 

374.000 

4.643 

.000 

4.394 

.110 

95.993 

.648 

376 

375.000 

4.652 

.000 

4.405 

.110 

96.242 

.650 

377 

376.000 

4.661 

.000 

4.415 

.110 

96.488 

.652 

378 

377.000 

4.670 

.000 

4.425 

.109 

96.734 

.653 

379 

378.000 

4.678 

.000 

4.436 

.109 

96.977 

.655 

380 

379.000 

4.687 

.000 

4.446 

.109 

97.219 

.657 

381 

380.000 

4.696 

.000 

4.456 

.109 

97.459 

.658 

382 

381.000 

4.704 

.000 

4.466 

.109 

97.698 

.660 

383 

382.000 

4.713 

.000 

4.476 

.108 

97.936 

.661 

384 

383.000 

4.722 

.000 

4.486 

.108 

98.172 

.663 

385 

384.000 

4.730 

.000 

4.496 

.  108 

98.406 

.665 

386 

385.000 

4.739 

.000 

4.507 

.108 

98.640 

.666 

387 

386.000 

4.748 

.000 

4.516 

.108 

98.872 

.668 

388 

387.000 

4.757 

.000 

4.526 

.107 

99.102 

.669 

389 

388.000 

4.765 

.000 

4.536 

.107 

99.332 

.671 

390 

389.000 

4.774 

.000 

4.546 

.107 

99.560 

.672 

391 

390.000 

4.783 

.000 

4.556 

.107 

99.787 

.674 

392 

391.000 

4.791 

.000 

4.566 

.106 

100.013 

.675 

393 

392.000 

4.800 

.000 

4.576 

.106 

100.238 

.677 

394 

393.000 

4.809 

.000 

4.585 

.106 

100.462 

.679 

395 

394.000 

4.817 

.000 

4.595 

.106 

100.685 

.680 

396 

395.000 

4.826 

.000 

4.605 

.106 

100.906 

.682 

397 

396.000 

4.835 

.000 

4.614 

.105 

101.127 

.683 

398 

397.000 

4.843 

.000 

4.624 

.105 

101.347 

.685 

399 

398.000 

4.852 

.000 

4.634 

.105 

101.566 

.686 

400 

399.000 

4.861 

.000 

4.643 

.105 

101.784 

.687 

401 

400.000 

4.870 

.000 

4.653 

.105 

102.001 

.689 

Figure  A5.  (Sheet  7  of  15) 


Appandix  A  Listing  of  Program  EQPR.FOR 


A35 


402 

401.000 

4.878 

.000 

4.662 

.  104 

102.218 

.690 

403 

402.000 

4.887 

.000 

4.672 

.  104 

102.434 

.692 

404 

403.000 

4.896 

.000 

4.681 

.104 

102 . 649 

.693 

405 

404.000 

4.904 

.000 

4.690 

.104 

102.863 

.695 

406 

405.000 

4.913 

.000 

4.700 

.104 

103 . 077 

.696 

407 

406.000 

4.922 

.000 

4.709 

.104 

103 . 290 

.698 

408 

407.000 

4.930 

.000 

4.718 

.103 

103 . 502 

.699 

409 

408.000 

4.939 

.000 

4.728 

.103 

103.714 

.700 

410 

409.000 

4.948 

.000 

4.737 

.103 

103.925 

.702 

411 

410.000 

4.957 

.000 

4.746 

.103 

104 . 136 

.703 

412 

411.000 

4.965 

.000 

4.755 

.  103 

104 . 347 

.705 

413 

412.000 

4.974 

.000 

4.764 

.  102 

104 . 557 

.706 

414 

413.000 

4.983 

.000 

4.773 

.102 

104 . 766 

.708 

415 

414.000 

4.991 

.000 

4.782 

.  102 

104.975 

.709 

416 

415.000 

5.000 

.000 

4.791 

.102 

105.184 

.710 

417 

416.000 

5.008 

.000 

4.800 

.102 

105.392 

.712 

418 

417.000 

5.016 

.000 

4.809 

.101 

105.599 

.713 

419 

418.000 

5.024 

.000 

4.818 

.  101 

105.805 

.715 

420 

419.000 

5.032 

.000 

4.827 

.101 

106.011 

.716 

421 

420.000 

5.040 

.000 

4.836 

.101 

106.215 

.717 

422 

421.000 

5.048 

.000 

4.845 

.101 

106.419 

.719 

423 

422.000 

5.056 

.000 

4.854 

.101 

106.621 

.720 

424 

423.000 

5.064 

.000 

4.863 

.100 

106.823 

.721 

425 

424.000 

5.072 

.000 

4.871 

.100 

107 . 024 

.723 

426 

425.000 

5.080 

.000 

4.880 

.100 

107.225 

.724 

427 

426.000 

5.088 

.000 

4.889 

.100 

107.424 

.726 

428 

427.000 

5.096 

.000 

4.897 

.100 

107.623 

.727 

429 

428.000 

5.104 

.000 

4.906 

.099 

107.821 

.728 

430 

429.000 

5.112 

.000 

4.915 

.099 

108.019 

.730 

431 

430.000 

5.120 

.000 

4.923 

.099 

108.216 

.731 

432 

431.000 

5.128 

.000 

4.932 

.099 

108.413 

.732 

433 

432.000 

5.136 

.000 

4.940 

.099 

108.606 

.734 

434 

433.000 

5.144 

.000 

4.949 

.099 

108.804 

.735 

435 

434.000 

5.152 

.000 

4.957 

.098 

108.999 

.736 

436 

435.000 

5.160 

.000 

4.966 

.098 

109.193 

.737 

437 

436.000 

5.168 

.000 

4.974 

.098 

109.387 

.739 

438 

437.000 

5.176 

.000 

4.983 

.098 

109.581 

.740 

439 

438.000 

5.184 

.000 

4.991 

.098 

109.774 

.741 

440 

439.000 

5.192 

.000 

4.999 

.097 

109.967 

.743 

441 

440.000 

5.200 

.000 

5.008 

.097 

110.159 

.744 

442 

441.000 

5.208 

.000 

5.016 

.097 

110.352 

.745 

443 

442.000 

5.216 

.000 

5.024 

.097 

110.544 

.747 

444 

443. OCO 

5.224 

.000 

5.032 

.  097 

110.735 

.748 

445 

444.000 

5.232 

.000 

5.041 

.097 

110.927 

.749 

446 

445.000 

5.240 

.000 

5.049 

.096 

111.118 

.750 

447 

446.000 

5.248 

.000 

5.057 

.096 

111.309 

.752 

44S 

447.000 

5.256 

.000 

5.065 

.  096 

111.  SCO 

.753 

449 

448.000 

5.264 

.000 

5.073 

.096 

111.691 

.754 

450 

449.000 

5.272 

.000 

5.081 

.096 

111.882 

.756 

451 

450.000 

5.280 

.000 

5.089 

.096 

112 . 072 

.757 

452 

451.000 

5.288 

.000 

5.097 

.095 

112 . 263 

.758 

453 

452.000 

5.296 

.000 

5.105 

.095 

112.454 

.760 

454 

453.000 

5.304 

.000 

5.113 

.095 

112 . 644 

.761 

455 

454.000 

5.312 

.000 

5.121 

.095 

112.325 

.762 

456 

455.000 

5.320 

.000 

5.129 

.095 

113.026 

.763 

457 

456.000 

5.328 

.000 

5.137 

.095 

113.217 

.765 

458 

457.000 

5.336 

.000 

5.145 

.094 

113.408 

.766 

459 

458.000 

5.344 

.000 

5.153 

.094 

113 . 599 

.767 

460 

459.000 

5.352 

.000 

5.161 

.094 

113 . 790 

.769 

461 

460.000 

5.360 

.000 

5.168 

.094 

113.981 

.770 

Figure  A5.  (Sheet  8  of  1 5) 


A36 


Appondix  A  Listing  of  Program  EQPR.FOR 


462 

461.000 

5.368 

.000 

5.176 

.094 

114.173 

.771 

463 

462.000 

5.376 

.000 

5.184 

.093 

114.365 

.772 

464 

463.000 

5.384 

.000 

5.192 

.093 

114.557 

.774 

465 

464.000 

5.392 

.000 

5.199 

.093 

114.750 

.775 

466 

465.000 

5.400 

.000 

5.207 

.093 

114.942 

.776 

467 

466.000 

5.408 

.000 

5.215 

.093 

115.135 

.778 

468 

467.000 

5.416 

.000 

5.222 

.093 

115.329 

.779 

469 

468.000 

5.424 

.000 

5.230 

.092 

115.523 

.780 

470 

469.000 

5.432 

.000 

5.237 

.092 

115.717 

.782 

471 

470.000 

5.440 

.000 

5.245 

.092 

115.912 

.783 

472 

471.000 

5.448 

.000 

5.253 

.092 

116.107 

.784 

473 

472.000 

5.456 

.000 

5.260 

.092 

116.303 

.786 

474 

473.000 

5.464 

.000 

5.268 

.092 

116.499 

.787 

475 

474.000 

5.472 

.000 

5.275 

.091 

116.696 

.788 

476 

475.000 

5.480 

.000 

5.282 

.091 

116.893 

.790 

477 

476.000 

5.488 

.000 

5.290 

.091 

117.091 

.791 

478 

477.000 

5.496 

.000 

5.297 

.091 

117.289 

.792 

479 

478.000 

5.504 

.000 

5.305 

.091 

117.489 

.794 

480 

479.000 

5.512 

.000 

5.312 

.090 

117.688 

.795 

481 

480.000 

5.520 

.000 

5.319 

.090 

117.889 

.796 

482 

481.000 

5.528 

.000 

5.326 

.090 

118.090 

.798 

483 

482.000 

5.536 

.000 

5.334 

.090 

118.292 

.799 

484 

483.000 

5.544 

.000 

5.341 

.090 

118.495 

.800 

485 

484.000 

5.552 

.000 

5.348 

.090 

118.698 

.802 

486 

485.000 

5.560 

.000 

5.355 

.089 

118.903 

.803 

487 

486.000 

5.568 

.000 

5.362 

.089 

119.108 

.804 

488 

487.000 

5.576 

.000 

5.370 

.089 

119.314 

.806 

489 

488.000 

5.584 

.000 

5.377 

.089 

119.521 

.807 

490 

489.000 

5.592 

.000 

5.384 

.089 

119.729 

.809 

491 

490.000 

5.600 

.000 

5.391 

.088 

119.938 

.810 

492 

491.000 

5.608 

.000 

5.398 

.088 

120.147 

.811 

493 

492.000 

5.616 

.000 

5.405 

.088 

120.358 

.813 

494 

493.000 

5.624 

.000 

5.412 

.088 

120.570 

.814 

495 

494.000 

5.632 

.000 

5.419 

.088 

120.783 

.816 

496 

495.000 

5.640 

.000 

5.426 

.087 

120.997 

.817 

497 

496.000 

5.648 

.000 

5.432 

.087 

121.212 

.819 

498 

497.000 

5.656 

.000 

5.439 

.087 

121.428 

.820 

499 

498.000 

5.664 

.000 

5.446 

.087 

121.645 

.822 

500 

499.000 

5.672 

.000 

5.453 

.087 

121.863 

.823 

501 

500.000 

5.680 

.000 

5.460 

.087 

122.083 

.825 

502 

501.000 

5.688 

.000 

5.467 

.086 

122.304 

.826 

503 

502.000 

5.696 

.000 

5.473 

.086 

122.526 

.828 

504 

503.000 

5.704 

.000 

5.480 

.086 

122.749 

.829 

505 

504.000 

5.712 

.000 

5.487 

.086 

122.974 

.831 

506 

505.000 

5.720 

.000 

5.493 

.085 

123.200 

.832 

507 

506.000 

5.728 

.000 

5.500 

.085 

123.427 

.834 

508 

507.000 

5.736 

.000 

5.507 

.085 

123.656 

.835 

509 

508.000 

5.744 

.000 

5.513 

.085 

123.886 

.837 

510 

509.000 

5.752 

.000 

5.520 

.085 

124.117 

.838 

511 

510.000 

5.760 

.000 

5.526 

.084 

124.350 

.840 

512 

511.000 

5.768 

.000 

5.533 

.084 

124.585 

.841 

513 

512.000 

5.776 

.000 

5.539 

.084 

124.821 

.843 

514 

513.000 

5.784 

.000 

5.546 

.084 

125.058 

.845 

515 

514.000 

5.792 

.000 

5.552 

.084 

125.297 

.846 

516 

515.000 

5.800 

.000 

5-559 

.083 

125.538 

.848 

517 

516.000 

5.808 

.000 

5.565 

.083 

125.780 

.850 

518 

517.000 

5.816 

.000 

5.571 

.083 

126.024 

.851 

519 

518.000 

5.824 

.000 

5.578 

.083 

126.270 

.853 

520 

519.000 

5.832 

.000 

5.584 

.082 

126.517 

.855 

521 

520.000 

5.840 

.000 

5.590 

.082 

126.766 

.856 

Figure  A5.  (Sheet  9  of  15) 


Appendix  A  Listing  of  Program  EQPR.FOR 


A37 


522 

521.000 

5.848 

.000 

5.596 

.082 

127 . 017 

.858 

523 

522.000 

5.856 

.000 

5.603 

.082 

127.269 

.860 

524 

523.000 

5.864 

.000 

5.609 

.082 

127.524 

.861 

525 

524.000 

5.872 

.000 

5.615 

.081 

127.780 

.863 

526 

525.000 

5.880 

.000 

5.621 

.081 

128.038 

.865 

527 

526.000 

5.888 

.000 

5.627 

.081 

128.298 

.867 

528 

527.000 

5.896 

.000 

5.633 

.081 

128.560 

.868 

529 

528.000 

5.904 

.000 

5.639 

.080 

128.823 

.870 

530 

529.000 

5.912 

.000 

5.645 

.080 

129.089 

.872 

531 

530.000 

5.920 

.000 

5.651 

.080 

129.357 

.874 

532 

531.000 

5.928 

.000 

5.657 

.080 

129.627 

.876 

533 

532.000 

5.936 

.000 

5.663 

.  079 

129.898 

.877 

534 

533.000 

5.944 

.000 

5.669 

.079 

130.172 

.879 

535 

534.000 

5.952 

.000 

5.675 

.079 

130.448 

.881 

536 

535.000 

5.960 

.000 

5.681 

.079 

130.726 

.883 

537 

536.000 

5.968 

.000 

5.687 

.078 

131. 007 

.885 

538 

537.000 

5.976 

.000 

5.692 

.078 

131.289 

.887 

539 

538.000 

5.984 

.000 

5.698 

.078 

131.574 

.889 

540 

539.000 

5.992 

.000 

5.704 

.078 

131.861 

.891 

541 

540.000 

6.000 

.000 

5.710 

.077 

132.150 

.893 

542 

541.000 

6.004 

.000 

5.715 

.077 

132.440 

.895 

543 

542.000 

6.009 

.000 

5.721 

.077 

132.729 

.896 

544 

543.000 

6.013 

.000 

5.726 

.077 

133 . 016 

.898 

545 

544.000 

6.017 

.000 

5.732 

.  076 

133.302 

.900 

546 

545.000 

6.022 

.000 

5.738 

.076 

133.587 

.902 

547 

546.000 

6.026 

.000 

5.743 

.076 

133.870 

.904 

548 

547.000 

6.030 

.000 

5.749 

.076 

134.153 

.906 

549 

548.000 

6.035 

.000 

5.754 

.075 

134.434 

.908 

550 

549.000 

6.039 

.000 

5.759 

.075 

134.714 

.910 

551 

550.000 

6.043 

.000 

5.765 

.075 

134.993 

.912 

552 

551.000 

6.048 

.000 

5.770 

.075 

135.271 

.914 

553 

552.000 

6.052 

.000 

5.776 

.074 

135.549 

.916 

554 

553.000 

6.057 

.000 

5.781 

.074 

135.825 

.917 

555 

554.000 

6.061 

.000 

5.786 

.074 

136.100 

.919 

556 

555.000 

6.065 

.000 

5.791 

.  074 

136.374 

.921 

557 

556.000 

6.070 

.000 

5.797 

.073 

136.647 

.923 

558 

557.000 

6.074 

.000 

5.802 

.073 

136.920 

.925 

559 

558.000 

6.078 

.000 

5.807 

.073 

137.191 

.927 

560 

559.000 

6.083 

.000 

5.812 

.073 

137.462 

.928 

561 

560.000 

6.087 

.000 

5.817 

.072 

137.732 

.930 

562 

561.000 

6.091 

.000 

5.823 

.072 

138.001 

.932 

563 

562.000 

6.096 

.000 

5.828 

.072 

138.270 

.934 

564 

563.000 

6.100 

.  000 

5.833 

.  072 

138.537 

.936 

565 

564.000 

6.104 

.000 

5.838 

.071 

138.804 

.937 

566 

565.000 

6.109 

.000 

5.843 

.071 

139 . 071 

.939 

567 

566.000 

6.113 

.000 

5.848 

.071 

139.336 

.941 

568 

567.000 

6.117 

.000 

5.853 

.071 

139 . 602 

.943 

569 

568.000 

6.122 

.000 

5.857 

.070 

139.866 

.945 

570 

569.000 

6.126 

.000 

5.862 

.070 

140.130 

.946 

571 

570.000 

6.130 

.000 

5.867 

.070 

140.393 

.948 

572 

571.000 

6.135 

.000 

5.872 

.070 

140.656 

.950 

573 

572.000 

6.139 

.000 

5.877 

.069 

140.919 

.952 

574 

573.000 

6.143 

.000 

5.882 

.069 

141.181 

.954 

575 

574.000 

6.148 

.000 

5.886 

.069 

141.442 

.955 

576 

575.000 

6.152 

.000 

5.891 

.  069 

141.704 

.957 

577 

576.000 

6.157 

.000 

5.896 

.069 

141.964 

.959 

578 

577.000 

6.161 

.000 

5.901 

.068 

142.225 

.961 

579 

578.000 

6.165 

.000 

5.905 

.068 

142.485 

.962 

580 

579.000 

6.170 

.000 

5.910 

.068 

142.745 

.964 

581 

580.000 

6.174 

.000 

5.915 

.068 

143.004 

.966 

Figure  A5.  (Sheet  10  of  15) 


A38 


Appendix  A  Listing  of  Program  EQPR.FOR 


582 

581.000 

6.178 

.000 

583 

582.000 

6.183 

.000 

584 

583.000 

6.187 

.000 

585 

584.000 

6.191 

.000 

586 

585.000 

6.196 

.  000 

587 

586.000 

6.200 

.000 

588 

587.000 

6.204 

.000 

589 

588.000 

6.209 

.000 

590 

589.000 

6.213 

.000 

591 

590.000 

6.217 

.000 

592 

591.000 

6.222 

.000 

593 

592.000 

6.226 

.000 

594 

593.000 

6.230 

.000 

595 

594.000 

6.235 

.000 

596 

595.000 

6.239 

.000 

597 

596.000 

6.243 

.000 

598 

597.000 

6.248 

.000 

599 

598.000 

6.252 

.000 

600 

599.000 

6.257 

.000 

601 

600.000 

6.261 

.000 

602 

601.000 

6.265 

.000 

603 

602.000 

6.270 

.000 

604 

603.000 

6.274 

.000 

605 

604.000 

6.278 

.000 

606 

605.000 

6.283 

.000 

607 

606.000 

6.287 

,000 

608 

607.000 

6.291 

.000 

609 

608.000 

6.296 

.000 

610 

609.000 

6,300 

,000 

611 

610.000 

6.304 

.000 

612 

611.000 

6.309 

.000 

613 

612.000 

6.313 

.000 

614 

613.000 

6,317 

.000 

615 

614.000 

6.322 

.000 

616 

615.000 

6.326 

.  000 

617 

616.000 

6,330 

.000 

618 

617.000 

6.335 

.000 

619 

618.000 

6,339 

.000 

620 

619.000 

6.343 

.000 

621 

620.000 

6.348 

.000 

622 

621.000 

6.352 

.000 

623 

622.000 

6.357 

.000 

624 

623.000 

6.361 

.000 

625 

624.000 

6.365 

.000 

626 

625.000 

6.370 

.000 

627 

626.000 

6.374 

.000 

628 

627.000 

6,378 

.000 

629 

628.000 

6,383 

.000 

630 

629.000 

6.387 

.000 

631 

630.000 

6.391 

.000 

632 

631.000 

6.396 

.000 

633 

632.000 

6.400 

.000 

634 

633.000 

6.404 

.000 

635 

634.000 

6.409 

.000 

636 

635.000 

6.413 

.000 

637 

636.000 

6.417 

.000 

638 

637.000 

6.422 

,000 

639 

638.000 

6.426 

.000 

640 

639.000 

6,430 

.000 

641 

640.000 

6.435 

.000 

5.919 

.067 

143.263 

.968 

5.924 

.067 

143.522 

.969 

5.928 

.067 

143.781 

.971 

5.933 

.067 

144.040 

.973 

5.937 

.066 

144.298 

.975 

5.942 

.066 

144.557 

.976 

5.946 

.066 

144.815 

.978 

5.951 

.066 

145.073 

.980 

5.955 

.065 

145.331 

.982 

5.959 

.065 

145.589 

.SOI 

5.964 

.065 

145.847 

.985 

5.968 

.065 

146.105 

.987 

5.972 

.065 

146.363 

.989 

5.977 

.064 

146.621 

.990 

5.981 

.064 

146.879 

.992 

5.985 

.064 

147.138 

.994 

5.989 

.064 

147.396 

.996 

5.994 

.063 

147.655 

.997 

5.998 

.063 

147.913 

.999 

6.002 

.063 

148.058 

1.000 

6.006 

.000 

.000 

.000 

.000 

.000 

.000 

.000 

.000 

.000 

.000 

.000 

.000 

.000 

.000 

.000 

.000 

.000 

.  000 

.000 

.000 

.000 

.000 

.000 

.000 

.000 

.000 

.000 

.000 

.000 

.000 

.000 

.000 

.000 

.000 

.000 

.000 

.000 

.000 

.000 

,000 

.000 

.000 

,000 

.000 

.000 

.000 

,000 

.000 

.000 

.000 

.000 

.000 

.000 

.000 

.000 

.000 

,  000 

,  000 

.000 

.000 

.000 

.000 

.000 

.000 

.000 

.000 

.000 

,000 

.000 

.000 

.000 

.000 

.000 

.000 

.000 

.000 

.000 

.000 

.000 

.000 

.000 

.  000 

.000 

.000 

.000 

.000 

.000 

.000 

.000 

.000 

.000 

.000 

.000 

.  000 

.000 

.000 

.000 

.000 

.000 

.000 

.000 

.000 

.000 

.000 

.000 

.000 

.000 

.000 

.000 

.  000 

.000 

.000 

.000 

.000 

.000 

.000 

.000 

.  000 

.000 

.000 

.000 

.000 

.000 

.000 

.000 

.000 

.000 

.000 

.000 

.000 

.000 

.000 

.000 

.000 

,000 

.000 

.000 

.000 

.000 

.000 

.000 

.000 

.000 

.000 

.000 

.000 

.000 

.000 

.000 

.000 

.000 

.000 

.000 

.000 

.000 

,000 

.000 

.000 

.000 

Figure  A5.  (Sheet  11  of  15) 


Appendix  A  Listing  of  Program  EQPR.FOR 


A39 


642 

641.000 

6.439 

.  000 

.000 

.000 

.000 

643 

642.000 

6.443 

.  000 

.000 

.000 

.000 

644 

643.000 

6.448 

.000 

.000 

.000 

.  000 

645 

644.000 

6.452 

.000 

.000 

.000 

.  000 

646 

645.000 

6.457 

.000 

.000 

.000 

.  000 

647 

646.000 

6.461 

.000 

.000 

.000 

.000 

648 

647.000 

6.465 

.000 

.000 

.000 

.000 

649 

648.000 

6.470 

.000 

.000 

.000 

.000 

650 

649.000 

6.474 

.000 

.000 

.000 

.000 

651 

650.000 

6.478 

.000 

.000 

.000 

.000 

652 

651.000 

6.483 

.000 

.000 

.000 

.000 

653 

652.000 

6.487 

.000 

.000 

.000 

.000 

654 

653.000 

6.491 

.000 

.000 

.000 

.  000 

655 

654.000 

6.496 

.000 

.000 

.000 

.  000 

656 

655.000 

6.500 

.000 

.000 

.000 

.000 

657 

656.000 

6.504 

.000 

.000 

.000 

.000 

658 

657.000 

6.509 

.000 

.000 

.000 

.  000 

659 

658.000 

6.513 

.000 

.000 

.000 

.000 

660 

659.000 

6.517 

.000 

.000 

.000 

.  000 

661 

660.000 

6.522 

.000 

.000 

.000 

.  000 

662 

661.000 

6.526 

.000 

.000 

.000 

.000 

663 

662.000 

6.530 

.000 

.000 

.000 

.000 

664 

663.000 

6.535 

.000 

.000 

.000 

.  000 

665 

664.000 

6.539 

.000 

.000 

.000 

.  000 

666 

665.000 

6.543 

.  000 

.000 

.000 

.  000 

667 

666.000 

6.548 

.000 

.000 

.000 

.000 

668 

667.000 

6.552 

.000 

.000 

.000 

.  000 

669 

668.000 

6.557 

.000 

.000 

.000 

.000 

670 

669.000 

6.561 

.000 

.000 

.000 

.000 

671 

670.000 

6.565 

.000 

.000 

.000 

.000 

672 

671.000 

6.570 

.000 

.000 

.000 

.000 

673 

672.000 

6.574 

.000 

.000 

.000 

.000 

674 

673.000 

6.578 

.000 

.000 

.000 

.000 

675 

674.000 

6.583 

.000 

.000 

.000 

.000 

676 

675.000 

6.587 

.000 

.000 

.000 

.000 

677 

676.000 

6.591 

.000 

.000 

.000 

.  000 

678 

677.000 

6.596 

.000 

.000 

.000 

.  000 

679 

678.000 

6.600 

.000 

.ooc 

.000 

.000 

680 

679.000 

6.604 

.000 

.000 

.000 

.000 

681 

680.000 

6.609 

.000 

.000 

.000 

.  000 

682 

681.000 

6.613 

.000 

.000 

.000 

.  000 

683 

682.000 

6.617 

.000 

.000 

.000 

.000 

684 

683.000 

6.622 

.000 

.000 

.000 

.000 

685 

684.000 

6.626 

.000 

.000 

.000 

.  000 

686 

685.000 

6.630 

.000 

.000 

.000 

.000 

687 

686.000 

6.635 

.000 

.000 

.000 

.  000 

688 

687.000 

6.639 

.000 

.000 

.000 

.  000 

689 

688.000 

6.643 

.000 

.000 

.000 

.  000 

690 

689.000 

6.648 

.  000 

.000 

.000 

.  000 

691 

690.000 

6.652 

.000 

.000 

.000 

.000 

692 

691.000 

6.657 

.000 

.000 

.000 

.000 

693 

692.000 

6.661 

.000 

.000 

.000 

.000 

694 

693.000 

6.665 

.000 

.000 

.000 

,  000 

695 

694.000 

6.670 

.000 

000 

.000 

.000 

696 

695.000 

6.674 

.000 

.000 

.000 

.000 

697 

696.000 

6.678 

.000 

.000 

.000 

.000 

698 

697.000 

6.683 

.000 

..000 

.000 

.000 

699 

698.000 

6.687 

.000 

.000 

.000 

.  000 

700 

699.000 

6.691 

.000 

.000 

.000 

.000 

701 

700.000 

6.696 

.000 

.000 

.000 

.000 

Figure  A5.  (Sheet  12  of  15) 


A40 


.000 
.000 
.000 
.000 
.000 
.000 
.000 
.000 
.000 
.000 
.000 
.000 
.000 
.000 
.000 
.000 
.000 
.000 
.000 
.  000 
.000 
.000 
.000 
.000 
.000 
.000 
.000 
.000 
.000 
.000 
.000 
.000 
.000 
.000 
.000 
.000 
.000 
.000 
.000 
.000 
.000 
.000 
.000 
.000 
.000 
.000 
.000 
.000 
.000 
.000 
.000 
.000 
.000 
.000 
.000 
.000 
.000 
.000 
.000 
.000 


Appendix  A  Listing  of  Program  EQPR.FOR 


702 

701.000 

6.700 

.000 

.000 

.000 

.000 

-000 

703 

702.000 

6.704 

.000 

.000 

.000 

.000 

.000 

704 

703.000 

6.709 

.000 

.000 

.000 

.000 

.000 

705 

704.000 

6.713 

.000 

.000 

.000 

.000 

.000 

706 

705.000 

6.717 

.000 

.000 

.000 

.000 

.000 

707 

706.000 

6.722 

.000 

.000 

.000 

.000 

.000 

708 

707.000 

6.726 

.000 

.000 

.000 

.000 

.000 

709 

708.000 

6.730 

.000 

.000 

.000 

.000 

.000 

710 

709.000 

6.735 

.000 

.000 

.000 

.000 

.000 

711 

710.000 

6.739 

.000 

.000 

.000 

.000 

.000 

712 

711.000 

6.743 

.000 

.000 

.000 

.000 

.000 

713 

712.000 

6.748 

.000 

.000 

.000 

.000 

.000 

714 

713.000 

6.752 

.000 

.000 

.000 

.000 

.000 

715 

714.000 

6.757 

.000 

,000 

.000 

.  000 

.000 

716 

715.000 

6.761 

.000 

.000 

.000 

.  000 

.000 

717 

716.000 

6.765 

.000 

.000 

.000 

.000 

.000 

718 

717.000 

6.770 

.000 

.000 

.000 

.000 

.000 

719 

718.000 

6.774 

.000 

.000 

.000 

.000 

.000 

720 

719.000 

6.778 

.000 

.000 

.000 

.000 

.000 

721 

720.000 

6.783 

.000 

.000 

.000 

.000 

.000 

722 

721.000 

6.787 

.000 

.000 

.000 

.  000 

.000 

723 

722.000 

6.791 

.000 

.000 

.000 

.000 

.000 

724 

723.000 

6.796 

.000 

.000 

.000 

.000 

.000 

725 

724.000 

6.800 

.000 

.000 

.000 

.000 

.000 

726 

725.000 

6.804 

.000 

.000 

.000 

.000 

.000 

'•27 

726.000 

6.809 

.000 

.000 

.000 

.000 

.000 

728 

727.000 

6.813 

,000 

.000 

.000 

.000 

.000 

729 

728,000 

6.817 

.000 

.000 

.000 

.000 

.000 

730 

729.000 

6.822 

.000 

.000 

.000 

.000 

.000 

731 

730.000 

6.826 

.000 

.000 

.000 

,  000 

.000 

732 

731.000 

6.830 

.000 

.000 

.000 

.000 

.000 

733 

732.000 

6.835 

.000 

.000 

.000 

.  000 

.000 

734 

733.000 

6,839 

.000 

.000 

.000 

.000 

.000 

735 

734.000 

6.843 

.000 

.000 

.000 

.000 

.000 

736 

735.000 

6.848 

.000 

.000 

.000 

.  000 

.000 

737 

736.000 

6,852 

.000 

.000 

.000 

,  000 

.000 

738 

737.000 

6,857 

.000 

.000 

,000 

.000 

.000 

739 

738,000 

6.861 

.000 

.000 

.000 

.000 

.000 

740 

739.000 

6.865 

.000 

.000 

.000 

.000 

.000 

741 

740.000 

6.870 

.000 

.000 

.000 

.000 

.000 

742 

741.000 

6.874 

.000 

.000 

.000 

.000 

.000 

743 

742,000 

6.878 

.000 

.000 

.000 

.000 

.000 

744 

743.000 

6.883 

.000 

.000 

.000 

.  000 

.000 

745 

744.000 

6.887 

.000 

.000 

.000 

.  000 

.000 

746 

745.000 

6.891 

.000 

.000 

.000 

.000 

.000 

747 

746.000 

6.896 

.000 

.000 

.000 

.000 

.000 

748 

747.000 

6.900 

.000 

.000 

.000 

.000 

.000 

749 

748.000 

6.904 

.000 

.000 

.000 

.  000 

.000 

750 

749.000 

6.909 

.000 

.000 

.000 

.000 

.000 

751 

750.000 

6.913 

.000 

.000 

.000 

.000 

.000 

752 

751.000 

6.917 

.000 

.000 

.000 

.000 

.000 

753 

752.000 

6.922 

.000 

.000 

.000 

.  000 

.000 

754 

753.000 

6.926 

.000 

.000 

.000 

.  000 

.000 

755 

754.000 

6.930 

.000 

.000 

.000 

.  000 

.000 

756 

755,000 

6.935 

.000 

.000 

.000 

.  000 

.000 

757 

756.000 

6.939 

.000 

.000 

.000 

.000 

.000 

753 

757.000 

6,943 

.000 

.000 

.000 

.  000 

.000 

759 

758.000 

6.948 

.000 

.000 

.000 

.  000 

.000 

760 

759.000 

6.952 

.000 

.000 

.000 

.000 

.000 

761 

760.000 

6.957 

.000 

.ooc 

.000 

.000 

.000 

Figure  A5.  (Sheet  13  of  15) 


Appendix  A  Listing  of  Program  EQPR.FOR 


A41 


762 

761.000 

6.961 

.000 

.000 

.000 

.000 

.000 

763 

762.000 

6.965 

.000 

.000 

.000 

.000 

.000 

764 

763.000 

6.970 

.000 

.000 

.000 

.000 

.000 

765 

764.000 

6.974 

.000 

.  000 

.000 

.000 

.  000 

766 

765.000 

6.978 

.000 

.000 

.000 

.000 

.000 

767 

766.000 

6.983 

.000 

.000 

.000 

.000 

.000 

768 

767.000 

6.987 

.000 

.000 

.000 

.000 

.000 

769 

768.000 

6.991 

.000 

.000 

.000 

.000 

.000 

770 

769.000 

6.996 

.000 

.000 

.000 

.000 

.000 

771 

770.000 

7.000 

.000 

.000 

.000 

.000 

.000 

772 

771.000 

7.004 

.000 

.000 

.000 

.000 

.000 

773 

772.000 

7.009 

.000 

.000 

.000 

.000 

.000 

774 

773.000 

7.013 

.000 

.  000 

.000 

.000 

.  000 

775 

774.000 

7.017 

.000 

.000 

.000 

.000 

.000 

776 

775.000 

7.022 

.000 

.000 

.000 

.000 

.000 

777 

776.000 

7.026 

.000 

.000 

.000 

.000 

.000 

778 

777.000 

7.030 

.000 

.000 

.000 

.000 

.000 

779 

778.000 

7.035 

.000 

.000 

.000 

.000 

.000 

780 

779.000 

7.039 

.000 

.000 

.000 

.000 

.000 

I 

DYEQ(I) 

VGENV(I) 

VUSED{1) 

PSV(I) 

1 

100.000 

324.335 

28.362 

.541 

2 

105.000 

299.957 

29.699 

.500 

3 

110.000 

276.084 

32.269 

.460 

4 

115.000 

252.851 

36.341 

.421 

5 

120.000 

230.273 

42.157 

.384 

6 

125.000 

208.366 

49.950 

.347 

7 

130.000 

187.145 

60.465 

.312 

8 

135.000 

166.631 

75.781 

.278 

9 

140.000 

146.844 

182.815 

.245 

10 

137.500 

156.686 

87.965 

.261 

11 

138.750 

151.778 

101.170 

.253 

12 

140.000 

146.844 

182.830 

.245 

13 

139.375 

149.307 

124.819 

.249 

14 

139.688 

148.112 

146.663 

.247 

15 

140.000 

146.844 

182.835 

.245 

16 

139.844 

147.511 

163.600 

.246 

17 

139.688 

148.112 

146.667 

.247 

18 

139.766 

147.812 

154.885 

.246 

19 

139.727 

147.962 

150.690 

.247 

20 

139.688 

148.112 

146.666 

.247 

21 

139.707 

148.037 

148.642 

.247 

22 

139.697 

148.075 

147.633 

.247 
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Detailed  Description  of  Program 
EQPR.FOR  and  Input  and 
Output  Files 


Introduction 


The  FORTRAN  program  EQPR.FOR  accepts  one  input  file  (EQPR.INP) 
specifying  the  characteristics  of  the  original  profile  and  the  essential 
characteristics  of  the  nourishment  material  and  placement  geometry.  The 
program  also  generates  one  output  file;  EQPR.OUT.  The  program  and  the 
input  and  output  files  are  described  in  this  appendix.  A  program  listing 
(EQPR.FOR),  and  annotated  input  and  output  files  for  Examples  1  and  4  have 
been  presented  in  Appendix  A.  In  the  following  presentation,  it  is  necessary 
to  recall  that  depths  and  elevations  represent  positive  and  negative  values  of 
depth  h,  respectively. 

Description  of  FORTRAN  Program  and  Input  Files 


This  discussion  of  the  FORTRAN  program  will  be  assisted  by  referring  to 
Part  2  and  Figure  10  of  the  main  body  of  this  report,  the  latter  of  which  is 
also  included  in  this  appendix  as  Figure  B1  for  reference. 

Main  program  EQPR.FOR 

The  primary  function  of  the  main  program  is  to  call  the  various  subroutines 
that  will  be  described  in  subsequent  sections  of  this  appendix.  In  addition,  in 
the  iterations  phase  of  this  program,  improved  estimates  of  the  equilibrium 
berm  displacement,  DYEQ,  are  generated.  The  overall  sequence  of  the  main 
program  is  as  follows. 

a.  Read  and  write  a  run  description  in  FORMAT  (20A4) 
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b.  Call  Subroutine  PRCXZHAR  which  reads  in  the  initial  and  placed 
profile  characteristics.  Solve,  by  iteration,  for  Ay„,  the  initial  berm 
displacement.  The  initial  profile  can  be  idealized  (of  the  following 
form) 


3 


or  may  be  user-specified. 

c.  Call  Subroutine  SEDCHAR,  which  reads  in  the  sediment  size 
information  in  the  form  of  an  idealized  (log-normal)  distribution  or  a 
user-specified  distribution. 

d.  Call  Subroutine  ATRANS,  which  converts  the  sediment  size 
distribution  information  to  a  distribution  of  the  sediment  scale 
parameter,  A.  This  program  reads  in  the  desired  number  of  iterations, 
ITSMAX,  and  an  initial  estimate  of  the  equilibrium  profile 
displacement,  DYEQl. 

e.  Call  Subroutine  PSSR,  which  manages  the  iteration  process,  as  will  be 
described  when  discussing  this  subroutine. 

g.  The  main  program  develops  improved  estimates  of  equilibrium  profile 
displacement  DYEQ. 

h.  Finally,  the  main  program  generates  output  to  FILE  EQPR.OUT, 
which  summarizes  the  results  of  the  run. 

The  paragraphs  below  describe  each  of  the  subroutines. 

Subroutine  PROCHAR.  This  subroutine  reads  in  the  characteristics  of  the 
initial  and  placed  profiles.  The  first  line  of  input  includes  the  following 
variables  in  FORMAT  (2I6,8F8.3);  JPTYPE,  IMAX,  AO,  SO,  SP,  VADD, 
SEQ,  HBERM,  HSTAR,  HMIX. 

JPTYPE  Type  of  profile  input.  If  JPTYPE  =  0,  the  depths  are  calculated  as 
an  equilibrium  profile  using  Equation  4.  If  JPTYPE  =  1,  the  depths  are  read 
in  as  individual  points  at  arbitrary  spacings  across  the  profile. 

IMAX  =  Number  of  final  points  across  the  profile. 

AO  =  Profile  scale  parameter  of  initial  profile  (only  required  if 

JPTYPE  =  0). 

DY  =  Spacing  of  points  across  final  profile. 

SO  =  Beach  face  slope  of  initial  parameter  (only  required  if  JPTYPE 

=  0). 

SP  =  Slope  of  placed  material. 

VADD  =  Volume  of  added  material  in  cubic  meters/meter. 
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SEQ  = 
HBERM  = 
HSTAR  = 

HMIX  = 


Slope  of  above  water  portion  of  equilibrium  profile. 

Height  of  berm  in  meters. 

Depth  of  effective  motion,  i.e.,  the  depth  to  which  the 
nourishment  material  is  considered  to  equilibrate. 

The  depth  to  which  nourished  sand  is  mixed  in  the  placement 
area. 


This  completes  information  contained  on  the  first  card.  The  program  then 
calls  Subroutine  PROCAL  if  JPTYPE  =  0  to  calculate  the  initial  profile.  If 
JPTYPE  =  1 ,  the  program  reads  in  IMAXP,  the  number  of  pairs  of  points  of 
(YV.HVO)  to  be  provided  as  input.  These  values  are  read  in  FORMAT 
(8F8.2).  Appendix  A  has  provided  annotated  input  files  for  Examples  1  and 
4,  respectively,  which  are  for  idealized  and  user-specified  initial  profiles.  The 
user-specified  profile  (YV,HV0)  is  then  interpolated  by  Subroutine  INTERP  to 
establish  IMAX  pairs  of  (Y,H0)  at  spacing  DY. 

With  the  information  above  describing  the  initial  profile  and  nourishment 
characteristics,  the  initial  shoreline  advancement  due  to  nourishment  is 
calculated.  The  characteristics  of  this  profile  are  a  berm  at  height  HBERM 
extending  a  distance  DYO  from  the  initial  shoreline,  sloping  at  SP  to 
intersection  with  the  initial  profile  (HO)  and  containing  VADD  volume  of 
sediment.  The  shoreline  advancement  DYO  is  determined  by  iteration. 

This  completes  description  of  the  Subroutine  PROCHAR. 

Subroutine  SEDCHAR.  This  subroutine  reads  in  the  characteristics  of  the 
nourishment  material.  The  first  card  has  the  following  format  (I6,2F8.2): 


JSED  = 


LMAX  = 


XMU  = 
SIG  = 


Type  of  sediment  input.  If  JSED  =  0,  the  sediment 
characteristics  are  described  by  Equation  1  and  XMU,  the 
sand  size  mean  in  phi  units  and  SIG,  the  sample  standard 
deviation,  also  in  phi  units.  If  JSED  =  1,  the  sediment 
characteristics  are  read  in  as  individual  pairs  of  (diameter, 
cumulative  distribution). 

Number  of  pairs  of  (sand  size,  decimal  percentage)  to  be 
generated  (JSED  =  0)  or  read  in  (JSED  =  1)  to  define  the 
cumulative  sand  size  distribution. 

Sand  size  in  phi  units  (JSED  =  0). 

Sand  standard  deviation  in  phi  units  (JSED  =  0). 


The  second  through  sixth  cards  read  by  this  subroutine  establish  a  table  of 
19  pai:  if  (DI,AI)  values  which  represent  the  sediment  scale  parameter 
variation  AI  with  diameter  DI.  These  values  will  be  the  same  for  all  runs  and 
are  simply  a  means  of  inputting  the  empirical  information  in  Figure  B2. 


If  JSED  =  0,  the  cumulative  distribution  of  the  nourishment  sediment  is 
calculated  in  accordance  with  Equation  1.  If  JSED  =  1,  the  cumulative 
distribution  is  read  in  as  pairs  of  (D(L),  P(L),  L  =  1 ,  LMAX)  in  format 
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(8F8.2)  where  D(L)  is  the  diameter  in  millimeters  and  P(L)  is  the  probability 
in  decimal  percentage.  Note  P(l)  =  0,  and  /*(LMAX)  =  1 .0. 


The  reader  is  referred  to  Appendix  A  for  examples  of  input  to  SEDCHAR 
for  idealized  (Figure  A2)  and  user-specified  (Figure  A4)  sediment 
characteristics,  respectively.  This  complete  description  of  the  Subroutine 
SEDCHAR. 

Subroutine  ATRANS.  This  subroutine  simply  calculates  sediment  scale 
values  A  for  each  of  the  diameters  D  in  the  probability  distribution.  Thus  the 
distribution  of  P  =  f(D)  is  transformed  to  A  =  f(P),  where  P  is  the  decimal 
percentage  coarser  dian. 

At  this  stage,  the  initial  profile,  placed  profile,  and  sediment 
characteristics  (including  their  A  values)  are  established.  The  main  program 
then  reads  in  the  number  of  iterations  to  be  performed,  ITS  MAX,  and  an 
initial  estimate  for  the  equilibrium  shoreline  advancement,  DYEQ  (the 
estimate  is  termed  DYEQl). 

The  main  program  calls  Subroutine  PSSR,  which  carries  out  the  solution. 

Subroutine  PSSR.  This  subroutine  establishes  the  equilibrium  profile 
locally  consistent  with  the  nourishment  volume  mobilized  and  grain  size 
distribution.  In  this  process,  the  program  considers  three  regions  as  discussed 
below  and  as  illustrated  in  Figure  B2. 

Iterations  are  carried  out  until  the  volume  of  sand  mobilized  is  equal  to  that 
deposited  out  to  the  point  where  the  equilibrium  and  initial  profiles  intersect 
(Intersecting  Profiles)  or  the  equilibrium  profile  reaches  HSTAR  (Non¬ 
intersecting  Profiles).  Common  to  all  regions  is  that  the  equilibrium  profile  is 
advanced  seaward  by  increments  of  DY  and  at  each  advancement,  the  volume 
generated  from  the  placement  area,  VGEN,  is  calculated  (if  appropriate)  and 
the  volume  used  (VUSED)  is  also  calculated.  Also  in  advancing  seaward,  the 
equilibrium  profile  computations  consider  coarsest  sand  to  be  deposited  in  the 
more  landward  portions  of  the  profile  and,  as  the  profile  advances  still  farther 
seaward,  finer  and  finer  sediment  is  utilized. 

a.  Region  I  -  This  is  the  region  in  which  sediment  is  removed  from  the 
placement  volume.  Above  the  water,  the  slope  is  considered  uniform 
as  SEQ.  Throughout  this  region,  the  volume  generated  is  considered 
to  include,  in  addition  to  the  net  volume  removed  from  the  placed 
geometry,  a  depth  HMIX,  which  is  user-specified.  The  equilibrium 
profile  is  extended  to  intersection  with  the  placed  profile.  All  of  the 
material  removed  up  to  this  point  (including  the  contribution  of 
HMIX)  represents  VGEN.  The  only  volume  used  (VUSED)  up  to  this 
point  is  that  due  to  HMIX.  Seaward  of  this  intersection  point,  VGEN 
does  not  change  and  VUSED  increases  as  the  summation  of  the 
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elemental  contributions  between  the  equilibrium  profile  and  the 
underlying  profile,  which  could  be  the  placed  or  initial  profile. 

b.  Region  II  -  This  region  lies  between  the  intersection  points  of  the 
equilibrium  and  placed  profiles  and  the  initial  and  placed  profiles. 

c.  Regioa  III  -  This  region  extends  from  the  intersection  of  the  initial  and 
placed  profile  to  the  end  of  the  equilibrium  profile,  which,  as  noted 
previously,  can  occur  due  to  profile  intersection  or  the  equilibrium 
profile  reaching  the  depth  of  effective  motion. 

In  the  calculations  of  the  equilibrium  profile,  the  Subroutine  PSSR  calls  two 
Subroutines,  YNEW  and  VOL,  each  of  which  is  described  below. 

a.  Subroutine  YNEW  -  This  subroutine  advances  the  equilibrium  profile 
by  a  distance  DY.  In  carrying  out  the  calculations,  an  estimate  of  the 
current  A  value  is  necessary.  This  is  determined  from  an  estimate  of 
the  total  value  of  VGEN  for  the  profile  and  the  current  value  of 
VUSED.  The  A  value  is  then  determined  from  A=f(P)  where  P  = 
VUSED/VGEN.  This  interpretation  is  accomplished  by  Subroutine 
INTERP. 

b.  Subroutine  VOL  -  This  subroutine  calculates,  depending  on  the  region 
of  equilibrium  profile  (vis-a-vis,  the  initial  and  placed  profiles),  the 
volume  increments  to  VGEN  and  VUSED  for  each  increment  of 
equilibrium  profile  advancement.  In  this  process,  it  is  extremely 
important  that  the  details  of  the  volumetric  contribution  be  evaluated 
correctly  at  the  intersection  points  (equilibrium  and  placed;  placed  and 
initial).  This  subroutine  calls  Subroutine  CROSS,  which  establishes 
the  distance  YINT  and  depth  HCROSS  of  the  intersection  point  of  the 
placed  and  equilibrium  profiles. 

The  output  file  EQPR.OUT  intbrmation  is  reasonably  well  defined  by  the 
output  annotation.  There  are  six  different  types  of  information  present^. 

The  first  line  is  simply  a  repeat  of  the  run  descriptor,  which  is  read  in  as  the 
first  card  in  EQPR.INP.  Next  are  the  input  pairs  of  Dl,  A1  also  provided  in 
the  input  files  to  all  runs.  Thirdly,  the  generated  or  input  pairs  of  sediment 
size  and  cumulative  probability  (D(I),  P(I))  are  presented.  The  next,  and  by 
far  the  largest  table,  includes  the  variation  of  the  following  with  distance,  Y; 
HOa),  HPG),  HEQO),  ACG),  VUSEDG),  and  PCCG),  where  all  terms  have 
been  described  previously,  except  for  PCCG),  which  is  the  value  of  the 
proportion  of  volume  used  up  to  that  point  relative  to  the  volume  available. 
Note  that  seaward  of  the  intersection  points  of  the  initial  (HOG))  and  placed 
(HPG))  profiles,  the  placed  profile  values  are  given  as  zero  and  similarly  for 
the  equilibrium  profile  seaward  of  intersection  with  the  initial  profile  or 
seaward  of  reaching  HSTAR.  The  last  table  presents  a  summary  of  the 
following  information  for  each  iteration  value:  DYEQG),  VGENG), 
VUSEDG),  and  PSVG).  Each  of  these  values  has  been  defined  before  with 
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the  exception  of  PSV(I),  which  is  the  ratio  of  total  volume  used  (and 
generated)  to  volume  added.  Note  that  for  the  solution  to  be  successful, 
VGEN  must  equal  VUSED.  The  last  line  of  the  output  file  announces  that  a 
solution  has  been  reached  and  whether  the  solution  is  of  the  intersecting 
profile  or  non-intersecting  profile  type. 

This  completes  the  detailed  description  of  the  computer  program 
EQPR.FOR,  the  input  file  EQPR.INP,  and  the  output  file  EQPR.OUT. 
Annotated  input  and  output  files  for  Examples  1  and  4  have  been  presented  in 
Appendix  A. 
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Figure  B1 .  Flow  diagram  for  problem  solution 
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Appendix  C 

Additional  Data  for  Delray 
Beach,  Florida 


This  appendix  presents  additional  beach  profile  and  grain  size  distribution 
data  for  Delray  Beach,  Florida.  Although  these  data  have  not  been 
analyzed/evaluated  to  date,  they  are  included  in  the  report  as  a  potentially 
valuable  data  source.  The  14  plots  are  described  briefly  below. 

a.  Figure  Cl  presents  the  average  grain  size  across  Delray  Beach  in 
1976.  These  results  are  for  a  composite  sample  based  on  several 
profiles.  Figure  C2  presents  the  grain  size  distributions  associated 
with  each  of  the  four  points  in  Figure  Cl. 

b.  Figures  C3,  C4,  C5,  and  C6  present  beach  profiles  for  pre¬ 
nourishment  (1973),  and  post-nourishment  (1983  and  1988)  for 
Stations  R177,  R180,  R184,  and  R187,  respectively. 

c.  Figures  C7,  C8,  C9,  and  CIO  present  grain  size  distributions  in  1988 
for  six  locations  across  the  profiles  at  Stations  R177,  R180,  R184,  and 
R187,  respectively. 

d.  Figures  Cll,  C12,  C13,  and  C14  present  the  1988  profiles  and 
associated  grain  size  distributions  at  three  locations  across  these 
profiles  at  Stations  R177,  R180,  R184,  and  R187,  respectively. 
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jure  C1 .  Average  grain  size  distribution  across  Delray  Beach,  Florida,  1976 
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Figure  C2.  Grain  size  distributions  at  various  locations  acros  the  profile,  Delray  Beach,  Florida,  1976 


Figure  C3.  Beach  profiles  at  R177  for  1973,  1983,  and  1988,  Di 
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igure  C6.  Beach  profiles  at  R187  for  1973,  1983,  and  1988,  Delray  Beach,  Florida 
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Grain  size  distributions  at  various  locations  across  profile  R177,  Delray  Beach,  Florida,  1988.  Grain  diameter  in 
millimeters 
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Figure  C8.  Grain  size  distributions  at  various  locations  across  profile  R180,  Delray  Beach,  Florida,  1988.  Grain  diameter  in 
millimeters 


C9.  Grain  size  distributions  at  various  locations  across  profile  R164,  Delray  Beach,  Florida,  1988.  Grain  diameter  in 
millimeters 
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Figure  C10.  Grain  size  distribution  at  various  locations  across  profile  R187,  Delray  Beach,  Florida,  1988.  Grain  diameter 
in  millimeters 
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Figure  Cl  2.  Beach  profile  at  Station  R180  in  1980  at  Delray  Beach,  Florida,  and  associated  grain  size  distributions.  Grain 
diameter  in  millimeters 
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Figure  Cl 4.  Beach  profile  at  Station  R187  in  1988  at  Delray  Beach,  Florida,  and  associated  grain  size  distributions.  Grain 
diameter  in  millimeters 
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A  Sediment  scale  parameter,  m*'^ 

Af  Nourishment  sediment  scale  parameter,  m*'^ 

Af,  Native  sediment  scale  parameter,  m‘'* 

B  Berm  height,  m 

D  Sediment  grain  size,  mm 

e  Euler’s  constant  (2.71828...) 

/  General  function;  grain  size  distribution 

F  Convergence  factor 

h  Water  depth,  m 

h„  Initial  profile  depth,  m 

/ip  Placed  profile  depth,  m 

h^Q  Equilibrium  profile  depth,  m 

h.  Depth  of  limiting  sediment  motion,  m 

Thickness  to  which  sediment  is  mixed  in  placed  portion  of 
cross  section,  m 

k  Iteration  number 

K  Empirical  coefficient,  m  ' 

m  Empirical  coefficient  describing  best-fit  variation  in  cross¬ 

shore  direction 

s  Beach  slope 

s„  Initial  profile  slope 

Sp  Nourishment  placement  slope 

V  Nourishment  volume  placed  per  unit  beach  length,  m^ 


D1 
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Sediment  fall  velocity,  m/sec 

Width  of  active  profile  for  the  initial  native  profile,  m 

Distance  seaward  tirom  the  shoreline 

Signifies  small  change  in  quantity 

Sediment  diameter  in  phi  units 

Sediment  sample  mean  in  phi  units 

Sediment  sample  deviation  in  phi  units 
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